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Introduction 


The  overall  goal  of  this  work  is  to  assess  the  potential  performance  characteristics  of  using  single  photon  emission 
computed  tomography  (SPECT)  for  imaging  inside  radiation  therapy  treatment  rooms  to  localize  functional  &  molecular 
targets  associated  with  metastatic  breast  cancer.  The  first  aim  is  to  evaluate  localization  accuracy  and  precision  as  a 
function  of  scan  time,  tumor  size,  tumor-to-background  uptake  ratio,  and  background  non-uniformity.  This  aim  was 
addressed  using  computer-simulation  and  hardware  studies.  The  second  aim  involves  the  localization  of  breast 
metastases  in  lung.  Through  computer  simulations,  we  studied  localization  performance  of  tumors  in  the  lungs  and 
other  sites  as  a  function  of  spatial  resolution  modeling,  observer  normalization,  iteration  number,  and  post¬ 
reconstruction  smoothing.  Localization  errors  were  also  analyzed  for  directional  dependencies  in  relation  to  the 
detector  trajectory.  Proposed  studies  on  respiratory  gating  and  the  use  of  compact  detectors  are  incomplete. 

Body 

Specific  Aim  1:  Evaluate  localization  accuracy  and  precision  as  a  function  of  scan  time,  tumor  size,  tumor-to- 
background  uptake  ratio,  and  background  non-uniformity. 

This  aim  was  addressed  by  assessing  localization  performance  using  (1)  computer  simulations  and  (2)  a  hardware 
study  with  a  whole-body  phantom  and  a  clinical  SPECT  scanner. 

(1)  Computer-simulation  study 


Fig.  1 .  Activity  distribution  (LEFT)  simulating  "mTc-sestarnibi  with  tumors  located  in  a  transaxial  slice  superior  to  the 
heart.  Three  tumor  diameters  and  two  uptake  ratios  were  investigated.  (RIGHT)  Attenuation  map  for  140  keV  photons 
from  "mTc  shown  with  anterior  half-circular  detector  orbit  [1], 

Tumors  were  simulated  with  various  diameters  -  10.8,  14.4  and  21.6  mm  -  and  with  uptake  ratios  of  3:1  and  6:1 
relative  to  background  and  were  located  in  bone,  lung  and  soft  tissue  of  a  female  NCAT  phantom  (Fig.  1).  These  tumor 
sizes  have  been  observed  clinically,  and  the  uptake  ratios  are  based  on  breast  cancer  imaging  with  Tc-99m  cardiac 
perfusion  radiotracers  [2], 

The  8  deepest  tumors  were  positioned  in  lung  and  were  subject  to  the  effects  of  respiratory  motion.  The  same  motion 
vectors  were  applied  to  laterally-sym metric  lung  tumors  in  order  to  minimize  differences  other  than  activity  uptake 
ratios.  Background  activity  was  non-uniform  for  tumors  near  the  heart.  An  anterior,  half-circular  detector  trajectory  was 
used  to  avoid  viewing  through  the  flat-top  treatment  couch.  The  simulated  projection  images  included  the  effects  of 
distance-dependent  spatial  resolution  and  non-uniform  attenuation.  Scatter  was  not  modeled.  Projections  were 
degraded  with  Poisson  Noise  to  form  ensembles  of  25  independent  noisy  realizations  for  scans  times  of  4,  8  and  20 
minutes.  Images  were  reconstructed  with  attenuation  correction  using  6-subset,  5-iteration  OSEM.  Representative 
reconstructed  images  are  shown  in  Fig.  2  from  the  image  ensembles. 


4 


dl 


-  f 


E 

ra  * 

TS  - 


o 

E 

3 


* 

« 

1  f 

% 

t 

to 

* 

»* 

<N 

1] 

VI 

i  % 

.  * 

scan  time  [min] 


i 

% 

* 

irr  • 


I  • 


t 

% 


*  1 


20 


4 

• 

% 

* 

* 

* 

• 

• 

§ 

• 

• 

% 

* 

i 

* 

• 

ft 

* 

> 

Fig.  2.  Sample  reconstructed  images  for  each  tumor  diameter  and  scan  time.  Grids  are  superimposed  on  images  such 
that  true  tumor  locations  are  centered  within  the  grid  squares  [1]. 


Tumors  were  localized  using  numerical  observers,  which  had  no  knowledge  of  noise  correlations,  for  a  multiple 
alternative  forced  choice  (MAFC)  task  where  one  and  only  one  tumor  was  located  within  each  28.8-mm-diameter 
search  volume.  Observer  templates  were  generated  by  blurring  true  tumor  geometries  with  a  14.4-mm-FWHM 
Gaussian  kernel,  which  approximates  image  spatial  resolution.  An  observer  was  translated  throughout  each  search 
volume  and  the  location  with  the  highest  cross  correlation  value  was  selected  as  the  measured  tumor  centroid. 
Euclidean  distances  between  the  measured  and  true  centroid  locations  were  computed  over  the  image  ensembles 
yielding  estimates  of  mean  error  --  localization  accuracy  -  and  standard  deviation  of  error  -  localization  precision. 
Estimates  of  localization  accuracy  and  precision  are  reported  in  Fig.  3  for  the  investigated  tumor  diameters,  anatomical 
locations,  tumor-to-background  activity  ratios  and  scan  times. 


Fig.  3.  Localization  accuracy  -  boxes  -  and  localization  precision  -  error  bars  -  displayed  as  a  function  of  tumor  size, 
anatomical  location,  tumor-to-background  activity  uptake  ratio  and  scan  time  [1], 
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The  results  in  Fig.  3  have  important  implications  for  using  on-board  SPECT  to  localize  metastatic  breast  cancer.  For 
stationary  tumors  of  the  same  diameter  and  uptake  ratio,  localization  accuracy  varied  substantially  depending  on  tumor 
location  relative  to  the  detector  trajectory.  Localization  accuracy  was  vastly  better  in  the  chest  wall  than  near  spine. 
These  differences  are  due  to  the  effects  of  attenuation  and  distance-dependent  spatial  resolution.  Because  of  the 
anterior  detector  trajectory,  relatively  fewer  counts  from  tumors  near  the  spine  escaped  the  phantom,  and  those  that 
did  escape  were  blurred  over  a  larger  number  of  detector  bins  as  compared  with  counts  from  tumors  near  the  chest 
wall.  While  the  investigated  half-circular  detector  trajectory  yields  highly  variable  localization  values  and  does  not 
appear  to  be  a  one-size-fits-all  acquisition  strategy,  it  is  important  to  note  that  the  task  of  on-board  SPECT  likely  would 
not  be  to  image  a  large  volume,  e.g.,  the  upper  torso.  Instead  the  task  would  be  to  image  a  relatively  small  region  near 
tumor.  Since  tumor  location  is  known  approximately  from  treatment  planning  images,  detector  trajectories  could  be 
customized  for  certain  anatomical  locations  to  minimize  the  effects  of  attenuation  and  distance-dependent  spatial 
resolution.  Further  research  is  needed  to  determine  impact  of  custom  detector  trajectories  on  target  localization,  given 
the  constraints  of  the  treatment  couch,  particularly  for  posterior  tumors  where  localization  was  generally  poor  using  an 
anterior  detector  trajectory. 

Results  from  Fig.  3  also  demonstrate  the  importance  of  radiotracer  specificity.  Tumors  in  laterally-symmetric 
anatomical  locations  differed  primarily  by  tumor-to-background  activity  ratios.  For  those  tumors  where  SPECT  was 
useful  for  localization,  tumors  of  6:1  uptake  were  typically  localized  as  well  using  4-minute  scans  as  tumors  of  3:1 
uptake  using  20-minute  scans.  Because  short  imaging  times  are  necessary  to  minimize  patient  motion  and  for 
machine  throughput,  radiotracers  with  a  high  differential  uptake  between  tumor  and  nearby  healthy  tissue  would  be 
most  appropriate  for  on-board  SPECT  imaging.  Additionally,  treatment  planning  SPECT  images  could  be  used  to 
estimate  tumor-to-background  activity  ratios  to  determine  which  patients  might  benefit  from  on-board  SPECT  guided 
radiation  therapy.  Planning  images  can  also  be  used  to  estimate  tumor  geometry,  another  important  criterion  for 
patient  selection  since  localization  accuracy  is  strongly  related  to  tumor  diameter,  as  shown  in  Fig  3. 

Non-uniform  background  activity  from  the  heart  proved  challenging  for  localizing  tumors  of  3:1  uptake  in  the  medial 
chest  wall.  These  tumors  are  denoted  by  black  boxes  in  the  first  column  of  Fig  3.  For  the  21 .6-mm-diameter  tumors, 
error  bars  are  much  smaller  than  the  total  error,  demonstrating  a  bias  that  can  be  visualized  in  reconstructed  images. 
Apparent  tumor  activity  is  suppressed  interiorly  near  the  heart  and  causes  the  measured  tumor  centroid  to  shift 
superiorly.  Additional  research  is  needed  to  improve  target  localization  near  hot  background  objects. 

Localization  was  generally  poor  in  the  lungs,  and  further  work  is  needed  to  investigate  the  impact  of  respiratory  gating 
on  localization.  As  expected,  localization  accuracy  and  precision  usually  improved  with  increasing  scan  time.  For  the 
shortest  scan  time,  certain  results  from  this  simulation  study  were  encouraging.  Localization  accuracy  was  within  2  mm 
for  anterior  tumors  of  6:1  uptake  that  were  14.4  mm  or  larger  in  diameter  using  4-minute  scans.  This  scan  time  is 
compatible  with  limited-fraction,  high-dose-per-fraction  therapies  like  stereotactic  body  radiation  therapy  where 
additional  time  is  allowed  for  patient  positioning. 

Further  details  pertaining  to  this  study  may  be  found  in  the  published  manuscript  [1],  which  is  located  in  Appendix  B. 
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(2)  Hardware  study 


In  related  SPECT  studies,  scatter  has  been  shown  to  degrade  observer  performance  in  the  detection  of  hot  tumors  in 
unknown  locations  [3].  The  above  computer-simulation  study  did  not  model  scatter  in  the  phantom  or  detector  [1].  In 
this  present  study,  localization  performance  was  assessed  in  scanner-acquired  SPECT  images  that  inherently  include 
scattered  photons  because  of  limited  energy  resolution.  This  study  is  important  for  gauging  the  degree  of  localization 
that  can  be  achieved  with  a  currently  available  SPECT  detector.  Also  investigated  in  this  study  are  the  effects  of  the 
detector  trajectory  on  localization  performance. 

The  hardware  phantom  consisted  of  five  fillable  spheres  -  targets  labeled  A-E  -  that  were  placed  in  a  single  40-cm 
section  of  the  Extended  Oval  PET  Phantom  [Data  Spectrum,  Hillsborough,  N.C.]  as  shown  in  Fig  4.  Inner  dimensions 
of  the  whole-body  phantom  were  34  cm  laterally  and  19  cm  anterior-posteriorly,  with  a  1-cm-thick  Acrylic  wall. 
Spherical  targets  were  weighed  before  and  after  filling,  and  by  assuming  a  density  of  1  g/cm3,  the  inner  diameters  of 
targets  A-E  were  calculated  as  2.8,  2.8,  3.4,  3.4,  and  2.2  cm.  Target-to-background  Tc-99m  radioactivity  ratio  was  6:1 
and  is  based  on  clinically  observed  tumor-to-background  activity  uptake  in  breast  cancer  using  cardiac  perfusion 
radiotracers  [2]. 

A 


E 

o 
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Fig.  4.  Shows  the  whole-body  phantom  and  five  fillable  spheres  labeled  A-E.  Inner  diameters  of  the  spheres  are 
reported  in  (cm). 

As  shown  in  Fig  5,  the  phantom  was  imaged  using  one  detector  of  a  Trionix  Triad  SPECT  scanner.  The  detector  was 
equipped  with  low  energy,  ultrahigh  resolution  (LEUR)  parallel-hole  collimation.  The  LEUR  collimator  holes  were  1.4 
mm  in  diameter  by  34.9  mm  in  length.  A  photopeak  window  for  "mTc  was  centered  on  140  keV  +/- 10.0%.  Projection 
images  were  acquired  every  2°  in  step-and-shoot  mode  and  implemented  on  a  256  x  128  grid  of  1 .78  x  1 .78  mm2  bins 


Fig.  5.  The  whole-body  phantom  was  imaged  by  one  detector  on  a  Trionix  Triad  SPECT  scanner  using  LEUR  parallel- 
hole  collimation. 
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Three  noncircular  detector  trajectories  were  investigated  -  left  lateral  180°,  360°,  and  right  lateral  180°  as  illustrated  in 
Fig  6  -  that  contour  to  the  phantom  and  couch,  thereby  achieving  better  spatial  resolution  than  is  possible  with  a  non- 
conformal  circular  detector  trajectory. 


Fig.  6.  Shows  the  three  conformal  detector  trajectories:  (a)  right  lateral  180°,  (b)  360°,  and  (c)  left  lateral  180° 

Scan  time  was  equivalent  for  these  detector  trajectories.  The  scan  time  was  adjusted  as  the  radioactivity  decayed  to 
maintain  the  equivalent  of  a  4.5-minute  clinical  scan  for  a  background  activity  level  of  0.25  pCi/ml.  In  regions  of  the 
detector  illuminated  by  background  activity,  approximately  7  photons  were  detected  in  the  photopeak  window  every 
second  per  cm2.  The  modeled  scan  time  does  not  include  the  time  spent  rotating  the  detector,  and  thus  models  the 
number  of  counts  expected  from  list-mode  acquisition,  which  could  be  implemented  for  on-board  SPECT  imaging. 

Ten  independent  scans  were  acquired  for  these  trajectories.  SPECT  images  were  reconstructed  on  a  grid  of  1 .78  x 
1.78  x  1.78  mm3  voxels  using  MLEM  without  corrections  for  attenuation,  distance-dependent  spatial  resolution,  or 
scatter. 

To  determine  localization  error,  the  true  position  of  a  target  must  be  known.  In  this  study,  the  true  positions  of  targets 
A-E  were  estimated  from  CT  and  12-hour  SPECT  images.  In  the  CT  image,  regions  of  interest  (ROIs)  were  placed  on 
the  spheres.  Intensity  thresholds  were  used  to  segment  air,  water,  and  sphere  wall.  Target  was  defined  as  sphere  wall 
and  internal  water,  excluding  air  pockets.  The  center  of  mass  was  calculated  for  each  target  and  was  used  as  an 
estimate  of  true  target  position  on  the  CT  grid.  Fig  7a  shows  the  CT  image  and  the  segmented  targets  averaged  over 
several  axial  slices.  In  a  12-hour  SPECT  image,  ROIs  were  placed  on  targets  A-E.  Background  activity  was  estimated 
from  proximate  slices  and  subtracted  from  each  ROI.  Target  appearance  in  SPECT  images  is  distorted  by  distance- 
dependent  spatial  resolution,  non-uniform  attenuation,  and  scatter.  Due  to  these  effects,  it  was  challenging  to  select  an 
optimal  intensity  threshold  based  on  visual  inspection  to  differentiate  target  from  background.  To  address  this 
challenge,  target  positions  in  SPECT  were  estimated  using  a  broad  range  of  intensity  thresholds.  Fig  7b  shows  the  12- 
hour  SPECT  image  and  the  targets  after  background  subtraction  and  intensity  thresholding. 


(b) 
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Fig.  7.  True  target  positions  were  derived  from  CT  and  12-hour  SPECT  images 
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The  CT-  and  SPECT-based  target  positions  -  5  pairs  of  3D  coordinates  -  were  rigidly  registered  using  Horn’s  closed- 
form,  least-squares  method  [4],  The  SPECT-based  position  estimates  with  the  lowest  residual  errors  were  considered 
the  truth  in  subsequent  assessments  of  localization  performance  in  the  4.5-minute  SPECT  images. 


Target  attenuation  and  distance  to  the  collimator  were  calculated  for  all  detector  views  in  each  detector  trajectory. 
Attenuation  was  estimated  from  the  CT  image  by  segmenting  air,  water,  and  Acrylic  and  by  assigning  140  keV  linear 
attenuation  coefficients  of  0.0,  0.155,  and  0.176  cm"r  respectively  to  these  materials.  Attenuation  survival  probability 
was  defined  as  the  fraction  of  photons  from  a  target  directed  towards  the  detector  that  were  not  attenuated  by  the 
phantom  or  couch.  Distance  was  calculated  from  a  target  centroid  to  the  patient-side  surface  of  the  collimator  along  a 
ray  perpendicular  to  the  detector  surface.  These  metrics  -  related  to  the  number  of  target  counts  and  collimator  spatial 
resolution  -  were  used  in  subsequent  comparisons  of  localization  performance  between  detector  trajectories. 

Target  localization  was  assessed  in  the  4.5-minute-equivalent  SPECT  images  using  a  numerical  observer.  Observer 
templates  were  binary  spheres  of  the  same  diameters  as  the  targets.  The  numerical  observer  searched  36-mm- 
diameter  volumes,  each  containing  only  one  target,  and  selected  the  location  with  the  highest  cross  correlation  value. 
Localization  error  was  the  distance  between  the  target  position  estimated  by  the  numerical  observer  and  the  target 
position  determined  from  the  12-hour  SPECT  and  CT  images. 


Fig.  8  shows  SPECT  images  from  the  three  detector  trajectories:  right  lateral  180°,  360°,  and  left  lateral  180°.  A  target 
is  most  readily  visible  in  an  image  corresponding  to  the  most  proximal  detector  trajectory,  i.e.,  targets  A  and  B  in  right 
lateral  180°  and  targets  D  and  E  in  left  lateral  180°.  Hot  focal  activity  near  the  phantom  surface  is  attributable  to 
external  fiducial  markers  not  used  in  this  analysis. 


Fig.  8.  Detector  trajectories  and  corresponding  SPECT  images:  (a)  right  lateral  180°,  (b)  360°,  and  (c)  left  lateral  180°. 
For  display  purposes,  images  were  summed  over  the  axial  slices  containing  at  least  one  target. 

In  Fig  9,  attenuation  survival  probability  and  distance  to  the  collimator  are  summarized  by  two  sets  of  curves,  each 
indexed  by  the  detector  angle.  Curves  on  the  top  row  show  the  distance  between  the  collimator  and  a  target.  Smaller 
values  indicate  better  spatial  resolution.  A  series  of  filled  curves  in  the  bottom  row  show  the  fraction  of  photons 
originating  from  the  center  of  a  target  that  are  not  attenuated  by  the  phantom  or  SPECT  couch.  The  area  under  these 
curves  is  proportional  to  number  of  tumor  counts  for  a  particular  detector  trajectory. 
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Fig.  9.  Collimator-target  distance  and  attenuation  survival  probability  are  plotted  as  function  of  the  detector  angle  and 
are  color-coded  by  detector  trajectory. 

The  distance  and  attenuation  results  displayed  in  Fig  9  were  averaged  over  each  detector  trajectory  and  are  reported 
in  Table  1 .  Target  C  was  centered  in  the  phantom  such  that  mean  attenuation  survival  probability  and  distance  to  the 
collimator  were  equivalent  for  each  trajectory.  For  the  other  targets,  a  proximal  180°  trajectory  yielded  the  lowest 
collimator-target  distances  and  the  highest  photon  survival  probabilities,  followed  by  360°  and  then  distal  180° 
trajectories. 

Table  1 :  Mean  target-collimator  distance  and  attenuation  survival  probability  for  each  detector  trajectory. 

Mean  Target-Collimator  Mean  Attenuation 

Distance  [cml_ Survival  Probability 


Target 

R  180° 

360° 

L  180° 

R  180° 

360° 

L  180° 

A 

11 

21 

30 

0.51 

0.31 

0.11 

B 

14 

21 

27 

0.34 

0.25 

0.17 

C 

21 

21 

21 

0.10 

0.10 

0.10 

D 

26 

21 

16 

0.09 

0.13 

0.18 

E 

29 

21 

13 

0.14 

0.26 

0.38 

Localization  errors  were  estimated  across  10  SPECT  images.  Results  are  displayed  in  Fig  10  using  box  plots.  Filled 
rectangles  include  data  within  the  25th  and  75th  percentiles.  Outliers  are  marked  by  filled  circles.  The  Wilcoxon  rank- 
sum  test  was  used  to  assess  for  statistical  significance  in  localization  differences  between  proximal  180°  and  360° 
trajectories.  The  symbol  *  indicates  p-value  <  0.05. 
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Fig.  10.  Ensemble  localization  errors  -  estimated  from  images  with  optimal  smoothing  and  iteration  number  -  are 
summarized  by  box  plots  and  color-coded  to  the  detector  trajectory.  The  symbol  *  indicates  that  localization  was 
significantly  better  using  a  proximal  180°  as  compared  with  a  360°  trajectory. 

Localization  errors  were  lowest  for  right-sided  targets  A  and  B  using  the  right  lateral  180°  trajectory,  were  slightly  larger 
with  the  360°  trajectory,  and  were  the  highest  and  substantially  so  with  the  left  lateral  180°  trajectory.  Similarly  for  left¬ 
sided  targets  D  and  E,  localization  errors  were  lowest  using  the  left  lateral  180°  trajectory  and  were  highest  for  the  right 
lateral  180°  trajectory.  These  results  demonstrate  that  localization  was  best  using  proximal  180°  trajectories. 
Localization  differences  between  proximal  180°  and  360°  trajectories  were  statistically  significant  at  A  and  D  with  p- 
values  of  0.004  and  0.049,  respectively.  Using  a  proximal  180°  trajectory,  mean  localization  error  was  within  1 .1  mm 
for  A,  B,  D,  and  E.  The  worst  localization  errors  across  the  10  images  were  within  2  mm  for  targets  A,  B,  and  E.  These 
results  demonstrate  that  a  commercially  available  SPECT  detector  can  provide  useful  information  for  localizing  certain 
targets  given  the  time  constraints  of  treatment  room  imaging. 

This  study  also  demonstrates  that  parallel-hole  SPECT  is  better  for  localizing  superficial  targets  than  deep  targets.  The 
volumes  of  deep  targets  C  and  D  were  1 .8  to  3.7  times  larger  than  superficial  targets  A,  B,  and  E;  however,  mean 
localization  errors  of  the  deep  targets  were  on  average  almost  twice  that  of  superficial  targets  when  using  a  proximal 
180°  trajectories.  Though  beyond  the  scope  of  this  study,  the  challenge  of  imaging  deep  targets  may  be  addressed  by 
using  converging  collimation. 
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Specific  Aim  2:  Evaluate  the  effect  of  respiratory  motion  on  the  accuracy  and  precision  of  target  localization. 

The  computer-simulation  study,  described  above  in  Aim  1,  modeled  respiratory  motion  for  tumors  in  the  lungs. 
Localization  performance  was  typically  poor  for  these  tumors,  as  shown  in  Fig  3.  Respiratory  motion  was  one 
degrading  effect  but  other  sources  of  poor  localization  performance  -  attenuation  and  spatially-varying  spatial 
resolution  -  are  inherent  to  SPECT  imaging.  In  light  of  findings  from  this  previous  study,  we  designed  another 
computer-simulation  study  to  analyze  localization  performance  for  lung  (and  other)  tumors  under  ideal  gating 
conditions,  i.e.,  no  motion,  and  for  a  scan  time  consistent  with  the  constraints  of  treatment  room  imaging. 

This  present  computer-simulation  study  expanded  upon  the  previous  computer-simulation  study  by  assessing  for 
directional  and  locational  dependent  localization  errors.  Also  investigated  were  the  effects  of  observer  normalization, 
detector  response  function  compensation,  iteration  number  and  post-reconstruction  smoothing  on  localization 
performance.  This  study  has  been  accepted  for  publication  in  the  peer-reviewed  journal  Medical  Physics.  For 
additional  details  on  this  study,  refer  to  the  corresponding  manuscript  in  Appendix  C. 

A  female  XCAT  Phantom,  an  extended  version  of  the  NCAT  Phantom  [5],  was  implemented  on  a  grid  of  0.15  x  0.15  x 
0.15  cm3  voxels.  Twelve,  1.5-cm-diameter  tumors  were  simulated  -  in  axilla,  breast,  chest  wall,  lung,  rib,  vertebra,  and 
other  soft  tissues  -  for  a  broad  sampling  of  locations  deep,  superficial,  anterior,  and  posterior  in  an  axial  slice  superior 
to  the  heart.  Relative  activity  ratios  were  20:6:1  for  heart,  tumor,  and  other  tissues  [2],  The  simulated  activity 
distribution  and  the  12  tumor  sites  labeled  A-L  are  shown  in  Fig  11a,b.  Respiratory  and  cardiac  motions  were  not 
modeled.  A  conformal  detector  trajectory  was  defined  such  that  for  each  detector  view,  the  detector  radius  was 
adjusted  to  achieve  0.5-cm  proximity  between  the  patient-side  surface  of  the  collimator  and  the  phantom  or  couch.  The 
detector  trajectory  is  illustrated  in  Fig  11c. 
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Fig.  11.  (a)  coronal  MIP  of  the  activity  phantom,  (b)  axial  slice  containing  12 
detector  trajectory  overlaid  on  the  attenuation  map. 


hot,  1 .5-cm-diameter  tumors,  and  (c)  the 


Images  were  reconstructed  using  OSEM  [6],  with  and  without  detector  response  function  compensation  (DRC)  [7],  and 
were  then  smoothed  with  a  3D  Gaussian  filter.  For  each  tumor  location  and  each  combination  of  DRC  versus  nDRC 
and  normalized  versus  un-normalized  observer  (defined  below),  localization  was  evaluated  across  a  broad  range  of 
image  updates  (i.e.  products  of  numbers  of  iterations  and  subsets)  and  of  smoothing  kernel  widths,  in  order  to 
determine  an  optimal  number  of  updates  and  degree  of  smoothing  at  each  location  and  for  each  method  of 
reconstruction  and  each  method  of  localization  estimation.  These  different  methods  were  compared  using  the  iteration 
numbers  and  smoothing  widths  that  were  optimal  for  each  method.  In  the  comparison  results  presented  below,  four 
subsets  were  used  to  accelerate  DRC  reconstructions,  but  only  one  subset  was  used  for  nDRC  images  because  the 
iteration-number  optimization  indicated  that  1  iteration  with  4  subsets  gave  too  many  image  updates  for  optimal 
localization  at  several  tumor  sites.  Post-reconstruction  smoothing  was  performed  using  3D  Gaussian  kernels  that 
ranged  in  FWHM  from  0  to  2.5  cm  in  increments  of  0.5  cm.  In  total,  80  noisy  images  were  reconstructed  using  2 
variations  of  OSEM  (DRC  and  nDRC)  and  smoothed  to  different  degrees  with  6  Gaussian  kernels  at  iterations  1-25. 


Localization  was  assessed  in  a  forced  choice  task,  where  each  3.0-cm-diameter  spherical  search  volume  Q.  contained 
one  and  only  one  tumor.  In  the  phantom,  the  tumor  was  centered  in  Q  at  coordinates  (x(,  yt,  zt).  A  numerical  observer 
sampled  Q  every  0.3  cm  in  Cartesian  space.  The  voxel  with  the  greatest  response  and  its  26  neighbors  were  then 
sampled  more  finely  at  every  0.1  cm.  The  observer  selected  subvoxel  location  i(xm  ,ym  ,zm)  within  Q  of  noisy  image  g 
that  produced  the  greatest  scalar  response  from  un-normalized  cross  correlation  (XC)  with  an  observer  template  w: 

z ( xm >  ym ’  zm)  =  arg max ( XC, ) ,  M) 

ieQ.  v  ' 
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where 


(2) 


XCi=wgr 

Thus,  location  (xm,  ym,  zm)  was  recorded  as  the  measured  tumor  position.  The  3D  template  w  was  a  sphere  -  the 
background-subtracted  tumor  signal  as  defined  in  the  XCAT  activity  phantom.  Template  dimensions  were  1 .5  x  1 .5  x 
1 .5  cm3.  The  same  template  w  was  used  for  all  1 2  tumor  sites  and  for  both  reconstruction  methods  and  variations  of 
cross  correlation.  Note  that  the  spherical  template  was  not  matched  exactly  to  a  tumor  because  tumor  appearance  in 
the  SPECT  images  was  distorted  due  to  the  effects  of  distance-dependent  spatial  resolution,  iteration  number,  and 
post-reconstruction  smoothing. 

Normalized  cross  correlation  (NXC)  was  also  investigated  for  estimating  tumor  positions: 

NXC;  =  ~j=  = '  =  (3) 

V(W-W)U«  •£») 

Localization  error  dr  was  calculated  as  the  Euclidean  distance  between  a  true  tumor  position  and  a  measured  tumor 
position: 

dr  =  ^t-xmf  +{.yt-ymf  +{zt-zmf  ■  (4) 

Localization  errors  were  recorded  across  the  noisy  image  ensembles  for  each  tumor  site.  Localization  errors  were 
compared  for  DRC-XC,  DRC-NXC,  nDRC-XC,  and  nDRC-NXC  using  site-specific  optimal  smoothing  and  iteration 
number  parameters. 


For  each  tumor  site,  direction-dependent  localization  errors  were  computed  every  3°  over  180°.  These  errors  were 
indexed  by  an  angle  6.  Angle  #was  defined  in  a  transaxial  plane  as  follows:  #was  equal  to  0  for  the  view  from  left 
lateral  to  center  of  rotation,  n/2  for  the  anterior  view,  and  n  for  the  right  lateral  view.  Direction-dependent  localization 

bias  (b^Q)  and  precision  indexed  to  a  specific  value  of  0,  were  calculated  along  the  dimension  perpendicular  to 
6.  As  such,  localization  precision  at  the  anterior  view  described  the  lateral  spread  of  localization  errors,  and 
localization  bias  -  the  lateral  bias.  Fig  12  illustrates  this  convention  using  data  from  tumor  G. 


The  direction-dependent  localization  errors  were  defined  mathematically  as 
follows:  Measured  tumor  coordinates  (xmk,  ymk )  from  noisy  image  k  and 

true  tumor  coordinates  (xt,  y, )  were  rotated  by  angle  fusing  matrix  T: 

cos 6  -sin# 

T=  (5' 

sin  9  cos  6 


to  yield 


and 


0 

X m,k 

=  T 

m,k 

_ym,k  _ 

ym,k 

e~\ 

r-  -| 

xt 

=  T 

ydt  _ 

Jt. 

(6) 

(7) 


For  the  direction  perpendicular  to  9,  direction-dependent  localization  bias 
Zf^was  calculated  as 


bw 


N 


N= 80 

X  iy^k  - 


k= 1 


(8) 


and  direction-dependent  localization  precision 


c^-0  as 


Fig.  12.  (a)  80  tumor  position  estimates 
are  marked  by  diamonds  and  overlaid 
on  the  attenuation  map.  The  dashed 
and  solid  circles  respectively  mark  the 
search  volume  and  tumor  boundaries. 

(b)  direction-dependent  localization  where 
precision  is  plotted  as  a  function  of 
angle  0. 
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Localization  performance  was  studied  across  a  broad  range  of  smoothing  filters  and  iteration  numbers.  Results  are 
reported  below  in  Fig  13.  Comparisons  in  localization  performance  between  nDRC-XC,  nDRC-NXC,  DRC-XC,  and 
DRC-NXC  were  made  using  the  combination  of  smoothing  and  iteration  number  that  resulted  in  the  lowest  mean 
localization  error.  Localization  errors  from  the  ensembles  of  80  images  are  summarized  by  box  plots  in  Fig  14. 


Localization  as  Function  of 
Smoothing  &  Iteration  Number 


Iteration  Number 


F° 

Tumor  Location  in  Phantom 


Fig.  13.  Shows  the  effects  of  smoothing  and  iteration  number  on  mean  localization  error  for  each  tumor  site  (A-L)  and 
also  for  the  variations  of  image  reconstruction  (DRC  vs.  nDRC)  and  observer  cross  correlation  (XC  vs.  NXC).  The 
combination  of  smoothing  and  iteration  number  that  yielded  the  lowest  mean  localization  error  is  marked  by  a  box. 
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Fig.  14.  Localization  errors,  calculated  across  ensembles  of  80  images,  are  summarized  by  box  plots.  Median  errors 
are  marked  by  horizontal  lines,  the  interquartile  range  (IQR)  by  filled  rectangles,  and  outliers  by  squares  or  diamonds. 
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The  random  chance  box  plot  describes  the  uniform  probability  of  selecting  any  location  in  a  search  volume  and  thus 
provides  a  benchmark  for  gauging  the  amount  of  information  provided  by  SPECT  for  target  localization. 

Tumor  position  estimates  are  displayed  on  the  attenuation  map  in  Fig  15.  Note  that  the  position  estimates  are  tightly 
clustered  at  site  B  but  are  more  diffuse  and  noticeably  anisotropic  at  site  G.  The  anisotropic  error  pattern  at  site  G  is 
related  to  the  detector  trajectory  in  Fig  16. 


Fig.  15.  Tumor  position  estimates  from  the  80-image  ensembles  are  overlaid  on  the  attenuation  map. 


Fig.  16.  (a)  Shows  the  detector  trajectory  in  reference  to  position  estimates  for  tumor  G.  (b)  Magnified  view  of  the 
position  estimates  demonstrates  that  detector  views  of  close  approach  and  minimal  attenuation  are  associated  with 
more  precise  localization. 

Direction-dependent  localization  bias  and  precision  -  calculated  using  Eqs  (8,9)  -  are  plotted  in  Fig  17  as  a  function  of 
angle  8.  Beneath  these  two  curves  are  two  additional  sets  of  curves:  The  first  set  of  curves  shows  collimator-tumor 
distance,  which  is  proportional  to  collimator  spatial  resolution,  and  the  second  set  shows  attenuation  survival 
probability  (defined  as  the  fraction  of  tumor  photons  directed  towards  the  detector  that  are  not  attenuated  by  the 
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phantom).  Detector  angles  of  close  approach  and  with  high  survival  probabilities  are  predictive  of  directions  with  the 
best  directional  bias  and  precision,  particularly  for  A-F. 
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Fig.  17.  (a)  Direction-dependent  localization  bias  and  precision  are  plotted  as  a  function  of  <9for  tumors  A-F.  Horizontal 
dashed  lines  mark  ideal  bias  and  precision.  Also  plotted  are  curves  showing  tumor-collimator  proximity  and  the  survival 
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probability  for  photons  from  a  tumor  to  escape  the  phantom  without  being  attenuated,  (b)  Shows  results  for  tumors  G- 
L.  Note  that  the  direction-dependent  localization  bias  curve  is  truncated  at  site  I. 

This  computer-simulation  study  compared  localization  performance  for  normalized  and  un-normalized  numerical 
observers,  which  were  used  to  estimate  tumor  positions  in  SPECT  images,  reconstructed  with  and  without  DRC.  For 
tumors  localized  to  <  2  mm  on  average,  which  are  good  candidates  for  SPECT-guided  radiation  therapy,  localization 
performance  typically  improved  by  compensating  for  the  detector  response  function  and  by  using  a  normalized 
observer. 

For  the  two  lung  tumors,  detector  response  function  compensation  did  not  improve  localization  performance.  In  certain 
images,  localization  errors  were  1  cm  or  greater  for  the  lung  tumors.  These  results  show  that  distance-dependent 
spatial  resolution  and  attenuation  can  have  a  major  effect  on  localization  performance.  Localization  likely  would  be 
worse  under  free-breathing  conditions.  Respiratory  gating  would  require  more  imaging  time  for  the  same  counting 
statistics.  These  considerations  and  the  results  above  demonstrate  the  challenges  to  imaging  lung  tumors  with 
parallel-hole  SPECT. 

Localization  errors  were  analyzed  across  the  image  ensemble  and  were  found  to  be  anisotropic  and  dependent  on 
tumor  position  relative  to  the  detector  trajectory.  Directional  dependencies  are  attributable  to  the  effects  of  attenuation 
and  distance-dependent  spatial  resolution.  Direction-dependent  localization  errors  have  important  implications  for 
external  beam  radiation  therapy  because  radiation  is  often  delivered  from  multiple  gantry  angles.  These  findings  on 
direction-dependent  localization  errors  might  be  useful  to  other  SPECT-guided  procedures,  e.g.,  breast  cancer  biopsy. 

The  proposed  studies  to  assess  the  localization  of  lung  tumors  as  a  function  of  different  respiratory  gating  strategies 
are  incomplete. 

Specific  Aim  3:  Evaluate  the  impact  of  limited  detector  surface  area,  and  consequent  projection-data  truncation, 
on  the  accuracy  and  precision  of  target  localization. 

Though  this  aim  is  largely  incomplete,  preliminary  results  in  Appendix  D,  demonstrate  that  a  sizeable  region  of  interest 
surrounding  a  tumor  can  be  reconstructed  with  similar  results  using  either  a  full-size  or  compact  detector  given  that  the 
conditions  of  [8]  are  met.  Localization  was  not  assessed  in  this  preliminary  study. 

Conclusions 

On-board  SPECT  was  investigated  for  localizing  functional  &  molecular  targets  associated  with  metastatic  breast 
cancer.  Computer-simulation  and  hardware  studies  were  performed  to  characterize  localization  accuracy  and  precision 
as  a  function  of  scan  time,  tumor  size,  tumor-to-background  uptake  ratio,  and  background  non-uniformity.  Localization 
performance  was  studied  in  images  with  and  without  spatial  resolution  modeling,  using  un-normalized  and  normalized 
observers,  and  for  a  broad  range  of  iteration  numbers  and  smoothing  parameters.  Of  the  parameters  that  can  be 
controlled,  the  detector  trajectory  typically  had  the  largest  effect  on  localization  performance.  A  proximal  180°  trajectory 
yielded  the  best  localization  performance  due  to  the  effects  of  distance-dependent  spatial  resolution  and  attenuation. 
Also  because  of  these  effects,  parallel-hole  SPECT  is  better  for  the  localization  of  superficial  targets  than  deep  targets. 
Another  parameter  that  can  be  controlled  is  observer  normalization.  For  tumors  proximal  to  hotter  background 
structures,  the  normalized  observer  placed  more  weight  on  template  shape  than  on  image  intensities,  as  compared 
with  the  un-normalized  observer,  and  in  these  cases,  localization  performance  was  vastly  better  using  a  normalized 
observer.  Mean  localization  errors  were  <  2  mm  for  superficial  tumors  of  >  1 .4  cm  in  diameter  and  of  6:1  uptake  ratio 
using  a  proximal  detector  trajectory  and  scan  time  of  approximately  5  minutes.  These  results  show  the  potential  of 
SPECT  for  localizing  certain  targets  given  the  time  constraints  of  imaging  inside  radiation  therapy  treatment  rooms. 
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Key  Research  Accomplishments 

Aim  1  was  addressed  through  computer  simulation  and  hardware  studies.  In  the  computer-simulation  study, 
localization  accuracy  and  precision  were  evaluated  as  a  function  of  scan  time,  tumor  size,  tumor-to-background  uptake 
ratio,  and  background  non-uniformity.  This  work  was  published  in  a  peer-reviewed  journal  Medical  Physics  (please  see 
Appendix  B.)  A  hardware  study  was  performed  to  study  the  effects  of  the  detector  trajectory  on  localization 
performance.  This  study  shows  the  degree  of  localization  accuracy  and  precision  that  can  be  achieved  using  a 
commercially  available  SPECT  detector.  Findings  from  this  study  have  been  submitted  to  Medical  Physics. 

Aim  2  involves  the  localization  of  tumors  in  the  lungs.  Localization  was  studied  for  these  and  other  tumors  as  a  function 
of  spatial  resolution  modeling,  observer  normalization,  number  of  iterations,  and  the  degree  of  post-reconstruction 
smoothing.  This  study  has  been  accepted  for  publication  in  Medical  Physics  (please  see  Appendix  C.) 

My  research  has  been  presented  at  conferences  hosted  by  AAPM  and  ASTRO.  Published  abstracts  are  listed  below  in 
Appendices  D-G. 

Other  Tasks  Related  to  BCRP  Predoctoral  Traineeship: 

(Appendix  A) 

•  I  have  gained  a  better  understanding  of  clinical  care  for  breast  cancer  patients  through  interactions  with  physicians 
from  radiology  and  radiation  oncology  at  Duke  University.  More  specifically,  I  have  attended  the  Joint  Breast 
Conference  where  patient  charts  are  reviewed  by  surgeons,  oncologists,  radiation  oncologists,  pathologists,  nurses, 
and  others  involved  in  patient  care  to  determine  the  best  course  of  treatment.  Further,  I  have  completed  the  Nuclear 
Medicine  Practicum  offered  by  the  Medical  Physics  Graduate  Program  where  I  spent  time  in  a  reading  room  with 
radiologists  who  explained  the  strategies  and  challenges  of  interpreting  nuclear  medicine  images. 

•  I  have  participated  in  and  presented  my  research  at  journal  clubs  and  seminars  hosted  by  the  Medical  Physics 
Graduate  Program  and  by  the  Department  of  Radiation  Oncology.  For  two  semesters,  I  was  responsible  for  organizing 
weekly  research  meetings  in  the  Physics  Division  of  Radiation  Oncology  where  ongoing  research  was  presented  by 
graduate  students  and  postdoctoral  researchers. 

•  I  have  completed  the  course  work  required  by  the  Medical  Physics  Graduate  Program  and  fulfilled  all  teaching 
assistantships.  During  the  spring  semester  2009,  I  participated  in  an  independent  study  involving  IMRT  treatment 
planning  under  the  supervision  of  a  senior  medical  physicist.  Through  this  course,  I  gained  practical  skills  for  creating 
treatment  plans  that  deliver  prescribed  dose  to  tumor  while  sparing  healthy  tissues. 

•  Results  from  this  report  were  included  in  my  preliminary  exam  proposal  and  presentation,  “On-board  SPECT  Imaging 
for  Localizing  Biological  Targets”,  which  was  accepted  by  my  committee  on  April,  20  2009. 

•  I  successfully  defended  my  dissertation  “On-board  Single  Photon  Emission  Computed  Tomography  (SPECT)  for 
Biological  Target  Localization”  on  June,  21  2010. 
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Appendix  A 


Specific  Aim  1:  Evaluate  localization  accuracy  and  precision  as  a  function  of  scan  time,  tumor  size,  tumor-to- 
background  uptake  ratio,  and  background  non-uniformity. 

A.  Hardware  Phantom  Studies  (Months  1-  8) 

a.  Fillable  hot  spheres  in  a  lower-activity  cylindrical  background,  approximately  mimicking  the 
geometry  of  supraclavicular  lymph  nodes 

b.  Sphere  diameters  ranging  from  8mm  to  24mm 

c.  Lesion-to-background  uptake  ratios  of  2:1,  5:1,  and  10:1 

d.  Background  non-uniformity:  40mm  sphere  with  an  offset  inner  8mm  sphere 

e.  Activity  ratios  of  10:2:1  and  10:5:1  in  the  outer  sphere,  inner  sphere,  and  background 

f.  Acquire  projection  images  using  a  Trionix  Triad  XLT  SPECT  scanner. 

g.  Scan  times  from  1  minute  to  24  minutes 

h.  Acquire  CT  image  of  the  same  phantom  (See  1 B  below). 

i.  From  ensemble  of  images,  evaluate  the  following  as  a  function  of  scan  time,  tumor  size,  tumor-to- 
background  uptake  ratio,  and  background  non-uniformity: 

i.  Evaluate  contrast-to-noise  ratios  (CNRs) 

ii.  Evaluate  MAFC 

j.  Prepare  and  submit  paper  on  this  work.  (Months  7-  8) 

B.  Use  Above  Hardware  Data  to  Validate  Computer-Simulation  Studies  (Months  9-12) 

a.  Generate  software  phantom  by  digitizing  the  CT  image  of  the  1A  hardware  phantom. 

b.  Repeat  1A  studies  using  the  digital,  software  phantom. 

c.  Compare  hardware  and  software  studies,  in  order  to  validate  the  software  approach  (which  is 
employed  below  in  2A-3A). 

Specific  Aim  2:  Evaluate  the  effect  of  respiratory  motion  on  the  accuracy  and  precision  of  target  localization. 

A.  Software  Studies  using  Digital  Anthropomorphic  Phantom  (Months  13-19) 

a.  Hot  spherical  lesions  surrounded  by  warm  background  activity  in  torso  of  NCAT 
anthropomorphic  phantom,  simulating  breast  cancer  metastasis  in  bone,  lung,  and  lymph  nodes. 

b.  Lesion  sizes  ranging  from  8mm  to  24mm 

c.  Radiotracer  uptake  ratios  of  2:1,  5:1,  and  10:1  in  lesions  and  background 

d.  Scan  times  ranging  from  1  minute  up  to  24  minutes 

e.  Generate  ensembles  of  breath-hold  images 

f.  Evaluate  CNR  and  MAFC  as  function  of  lesion  size,  lesion-to-background  uptake  ratio,  and  scan 
time. 

B.  Digital  Anthropomorphic  Phantom  with  Respiratory  Motion  (Months  20-  28) 

a.  Enable  respiratory  motion  features  within  NCAT  software. 

i.  Model  ungated  onboard  SPECT  imaging  during  free  breathing 

ii.  Model  gated  onboard  SPECT  imaging  during  free  breathing 

b.  Repeat  2A  studies  with  the  addition  of  respiratory  motion 

c.  Calculate  CNR  and  MAFC  as  a  function  of  lesion  size,  lesion-to-background  uptake  ratio,  scan  time, 
and  scan  mode  (as  described  in  2Bai-ii) 

d.  Prepare  and  submit  paper  for  work  in  2A-B.  (Months  27-  28) 

Specific  Aim  3:  Evaluate  the  impact  of  limited  detector  surface  area,  and  consequent  projection-data  truncation, 
on  the  accuracy  and  precision  of  target  localization. 

A.  Software  Studies  using  Truncated  Detector  Field  of  View  (FOV)  (Months  29-36) 

a.  Truncate  projection  images  from  2A  to  simulate  detectors  with  smaller  surface  areas. 

b.  Detector  FOV  ranging  from  15x15cm2  to  20x40cm2 

c.  Investigate  the  dependence  of  CNR  and  MAFC  on  detector  FOV. 

d.  Evaluate  the  relationship  between  detector  FOV  and  truncation  artifacts. 

e.  Write  and  submit  paper  (Months  35-36) 

Other  Tasks  Related  to  BCRP  Predoctoral  Traineeship. 

A.  Shadow  radiologist(s)  to  better  understand  breast  cancer  diagnoses  and  staging.  (Months  1-12) 

B.  Shadow  radiation  oncologist(s)  to  learn  about  breast  cancer  therapies.  (Months  1-12) 

C.  Participate  in  weekly  journal  clubs  and  seminars.  (Months  1-36) 

D.  Audit  courses  relevant  to  my  breast  cancer  research  project  (Months  1-36) 

E.  Present  research  at  conferences  such  as  AAPM,  ASTRO,  and  the  Era  of  Hope  Meeting.  (Months  1-36) 

F.  Gain  experience  with  RT  medical  physics  (e.g.  treatment  planning  &  quality  assurance).  (Months  12-36) 

G.  Prepare  and  defend  PhD  dissertation  (Months  30-36) 
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Single  photon  emission  computed  tomography  (SPECT)  was  investigated  for  imaging  on-board 
radiation  therapy  machines  in  order  to  localize  functional  and  molecular  targets.  A  computer- 
simulated  female  NCAT  phantom  was  positioned  supine  on  a  flat-top  treatment  couch.  Twenty  tumor 
locations  were  defined  in  the  upper  torso.  The  eight  lung  tumors  were  subject  to  the  effects  of 
respiratory  motion.  Tumor  diameters  of  10.8,  14.4,  and  21.6  mm  were  simulated  for  tumor-to- 
background  ratios  of  3:1  and  6:1  that  are  characteristic  of  the  radiotracer  "mTc-sestamibi.  Projec¬ 
tion  images  representing  scan  times  of  4,  8,  and  20  min  were  simulated  for  an  anterior,  half-circular 
trajectory.  Images  were  reconstructed  with  attenuation  correction  by  ordered- subsets  expectation 
maximization  (OSEM)  using  six  subsets  and  five  iterations.  Contrast-to-noise  ratios  (CNRs)  were 
calculated  from  ensembles  of  25  images.  Cross  correlation  with  a  noise-free  tumor  template  was 
used  to  select  the  most  suspicious  tumor  location  within  a  14.4-mm-radius  search  volume  surround¬ 
ing  each  tumor,  with  only  that  one  tumor  in  each  search  volume.  Localization  accuracy  was 
assessed  by  calculating  average  distances  between  measured  and  true  tumor  locations.  Localization 
accuracy  and  CNRs  were  strongly  affected  by  tumor  location  relative  to  the  detector  trajectory.  For 
example,  CNR  values  near  the  chest  wall  were  greater  by  a  factor  of  3.5  than  for  tumors  near  the 
spine  and  posterior  ribs,  a  much  greater  effect  than  the  factor  of  1.6  difference  in  CNR  between  6:1 
and  3:1  tumor  uptakes.  Typically,  tumors  of  6:1  uptake  were  localized  as  accurately  with  4  min  of 
scan  time  as  tumors  of  3:1  uptake  that  had  been  imaged  for  20  min.  Using  4  min  scans,  14.4  and 
21.6  mm  anterior  tumors  of  6:1  uptake  were  localized  within  2  mm.  These  results  suggest  that 
SPECT,  on-board  radiation  therapy  machines,  may  be  a  viable  modality  for  localizing  certain 
functional  and  molecular  targets  using  relatively  short  scan  times.  ©  2009  American  Association  of 
Physicists  in  Medicine.  [DOI:  10.1118/1.3113902] 

Key  words:  SPECT,  on-board  imaging,  target  localization 


I.  INTRODUCTION 

Functional  imaging  modalities  such  as  positron  emission  to¬ 
mography  (PET)  and  single  photon  emission  computed  to¬ 
mography  (SPECT)  rely  on  radiotracers  to  identify  physi¬ 
ological  processes — such  as  angiogenesis,  apoptosis, 
hormone  receptor  status,  hypoxia,  and  proliferation — that 
have  important  implications  in  cancer  management.1  In  the 
context  of  radiation  therapy,  such  functional  information  can 
be  used  to  define  biological  target  volumes  and,  when  com¬ 
bined  with  the  dose  sculpting  capabilities  of  intensity  modu¬ 
lated  radiation  therapy  (IMRT),  may  provide  enhanced  thera¬ 
peutic  success.2  To  date,  SPECT  and  PET  have  been  utilized 
for  diagnosis,  staging,  and  treatment  planning.  We  propose 
that  functional  imaging  may  also  have  an  important  role  in 
the  treatment  room  for  the  delivery  of  radiation  therapy.  Cur¬ 
rently,  treatment-room  imaging  is  composed  mainly  of  x-ray 
transmission  imaging  and  ultrasound  imaging,  which  are  pre¬ 
dominantly  anatomical  imaging  modalities.3  Imaging  with 
these  modalities  immediately  prior  to  delivery  has  improved 
patient  alignment  and  reduced  normal  tissue  damage  for 
some  sites,  but  there  are  limitations.  Notably,  anatomical 
differences  between  tumor  and  surrounding  tissue  may  not 
be  substantial  enough  to  directly  visualize  tumor.  Conse¬ 
quently,  localization  is  often  performed  instead  with  indirect 


markers  such  as  bony  anatomy,  that  may  not  necessarily  in¬ 
dicate  tumor  position.4-6  This  limitation  may  be  particularly 
important  when  high  gradients  are  used  to  sculpt  dose  around 
biological  targets  that  cannot  be  visualized  using  current 
treatment-room  imaging  modalities.  Further,  ultimately  ra¬ 
diation  therapy  may  be  directed  to  temporally  varying  func¬ 
tional  targets,  such  as  hypoxia,  which  at  treatment  time  may 
be  in  a  location  different  from  their  location  during  planning 
CT  or  planning  functional  imaging.7 

One  possible  method  for  directly  localizing  biological  tar¬ 
gets  is  to  image  with  SPECT  inside  radiation  therapy  treat¬ 
ment  rooms.  SPECT  has  been  used  for  a  wide  array  of  tumor 
imaging  applications,8,9  and  ongoing  research  is  expanding 
the  already  broad  array  of  radiotracers.10-15  While  spatial  res¬ 
olution  is  somewhat  limited  compared  with  many  anatomical 
imaging  modalities,  contrast  between  tumor  and  background 
can  be  substantially  greater,  and  this  stronger  tumor  signal 
may  improve  treatment-room  tumor  localization. 

Though  SPECT  has  been  studied  extensively  for  diagnos¬ 
tic  tasks,  SPECT  has  not  been  developed  for  the  treatment- 
room  environment  where  new  challenges  await.  Notably, 
SPECT  scans  typically  last  20  min  or  more,  and  this  amount 
of  time  is  too  long  in  terms  of  patient  motion  and  machine 
throughput.  Shorter  imaging  times  are  needed,  which  implies 
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Fig.  1.  (A)  Coronal  view  of  radiotracer  distribution  in  phantom;  (B)  uptake 
pattern  for  different  tumor  sizes  in  transaxial  slice. 


greater  image  noise.  Additionally,  treatment-room  couches 
are  wider  and  supporting  rails  are  typically  deeper  than  di¬ 
agnostic  SPECT  couches.  As  such,  treatment  couches  may 
challenge  detector-to-patient  proximity  and  increase  attenua¬ 
tion,  which  in  turn  degrade  image  quality  with  poorer  spatial 
resolution  and  increased  noise.  While  these  constraints  will 
negatively  impact  image  quality,  it  is  important  to  consider 
that  the  task  for  on-board  SPECT  would  be  different  than 
with  most  diagnostic  SPECT  scans.  Because  target  size  and 
approximate  location  are  known,  it  may  be  possible  to  accu¬ 
rately  localize  functional  targets.  Furthermore,  SPECT  detec¬ 
tor  trajectories  can  be  optimized  for  imaging  the  treatment 
region,  a  factor  which  may  enhance  SPECT  images  in  the 
region  of  interest  (ROI)  as  compared  to  diagnostic  scans 
which  typically  have  to  survey  a  larger  volume. 


II.  METHODS  AND  MATERIALS 
II.A.  Phantom 

A  female  anthropomorphic  phantom  was  generated  with 
NCAT  software.16,17  The  phantom  was  placed  supine  on  a 
flat-surface  treatment  couch.  Twenty  tumor  locations  were 
arranged  periodically  on  a  square  grid  with  centroid-to- 
centroid  spacing  of  approximately  6  cm.  Tumors  were  cen¬ 
tered  on  voxels  in  an  axial  slice  superior  to  the  heart.  For 
each  tumor  location,  three  tumor  diameters — 10.8,  14.4,  and 
21.6  mm  with  true  volumes  of  0.980,  2.10,  and  6.49  cm3 
— were  simulated  in  separate  phantoms  as  shown  in  Fig.  1. 
The  8  deepest  tumors  were  located  in  the  lungs  where  they 
were  subject  to  respiratory  motion,  while  the  remaining  12 
tumors  were  stationary  in  soft  tissues  or  bone.  Activity  dis¬ 
tribution  was  modeled  for  the  radiotracer  "mTc-sestamibi 
that  has  been  used  to  image  a  number  of  tumors,  including 
breast  cancer  and  its  metastases. 18-20  Activity  concentrations 
were  5.0  /zCi/ml  in  the  gall  bladder  and  myocardium, 
2.5  /zCi/ml  in  the  liver,  kidneys,  and  spleen,  and 
0.25  /zCi/ml  in  other  healthy  tissues.  Left-  and  right-sided 
tumor  activity  concentrations  were  1.75  and  1.5  /zCi/ml,  re¬ 


Fig.  2.  Attenuation  map  of  phantom  and  treatment  couch  with  dimensions. 


spectively,  as  these  uptake  values  relative  to  normal 
tissues — 3:1  and  6:1 — have  been  noted  clinically.19 

We  used  NCAT  software  to  generate  a  CT-like  image  of 
linear  attenuation  coefficients  for  140  keV  photons  in  the 
following  tissues:  0.149-0.157  cm-1  in  soft  tissues  and 
0.163-0.220  cm-1  in  bone.  The  attenuation  phantom  had  the 
exact  geometry  of  the  activity  phantom.  The  couch  was  as¬ 
signed  a  linear  coefficient  of  0.250  cm-1,  which  corresponds 
to  a  predominately  carbon-based  material  with  physical  den- 
sity  of  1.8  g/cnr.  A  map  of  attenuation  coefficients  and 
dimensions  for  a  slice  containing  tumor  are  shown  in  Fig.  2. 

The  effects  of  respiration  were  considered  by  averaging 
16  phases  of  a  respiratory  cycle  where  the  diaphragm  moved 
20.0  mm  and  the  chest  expanded  by  12.0  mm.  Normal  struc¬ 
tures  experienced  varying  degrees  of  2D  motion  in  the  an¬ 
teroposterior  (AP)  and  superoinferior  (SI)  dimensions.  Iden¬ 
tical  motion  vectors  were  applied  to  laterally  symmetric  left- 
and  right- sided  lung  tumors  to  minimize  differences  other 
than  activity  concentration.  Lung  tumor  motion  was  charac¬ 
terized  by  one  of  four  maximum  centroid-to-centroid  dis¬ 
placements  according  to  tumor  location:  anteromedial  (7.20 
mm  SI,  7.20  mm  AP),  anterolateral  (10.8  mm  SI,  10.8  mm 
AP),  posteromedial  (7.20  mm  SI,  3.60  mm  AP),  posterolat¬ 
eral  (10.8  mm  SI,  3.60  mm  AP).  Tumors  outside  the  lungs 
were  stationary.  Activity  and  attenuation  phantoms  were 
implemented  on  a  384X384X384  grid  of  0.18-cm-wide 
voxels.  Attenuation  phantoms,  one  for  each  tumor  size,  were 
also  placed  on  a  coarser  grid  with  0.36-cm-wide  voxels  and 
were  used  to  correct  for  nonuniform  attenuation  during  re¬ 
construction. 

II. B.  Projection  data  simulation  and  image 
reconstruction 

SPECT  imaging  was  simulated  with  an  analytical,  ray- 
driven  software  using  the  code  SPECT-MAP  which  has  been 
utilized  and  tested  in  this  and  many  previous  studies,  e.g., 
Refs.  22-26.  Projection  images  were  modeled  for  a  gamma 
camera  with  46X23  cm2  active  surface  area  and  parallel- 
hole,  low-energy,  high-resolution  collimation.  Simulated  col¬ 
limator  holes  were  2.7  cm  long  by  0.14  cm  in  diameter. 
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Fig.  3.  Line  source  profiles  for  measured  and  simulated  data  that  show  the 
effects  of  distance-dependent  collimator  resolution.  Measurements  were 
made  with  a  Trionix  Triad  SPECT  scanner. 


Distance-dependent  spatial  resolution  was  determined  by  fit¬ 
ting  an  analytical  form  to  measured  projection  data.  This 
spatial  resolution  modeling  was  implemented  by  tracing  a 
cone  of  121  rays  from  multiple  subdivisions  of  each  detector 
bin.  The  resolution  modeling  was  validated  by  comparing 
measured  line  spread  functions  to  those  calculated  using 
121 -ray  cones,  as  shown  in  Fig.  3.  Intrinsic  resolution  was 
modeled  with  a  2D,  0.34-cm-FWHM  Gaussian  kernel.  Also 
modeled  were  the  efficiencies  of  the  collimator,  scintillator, 
and  branching  ratio  of  "mTc.  SPECT-MAP  calculations  of  ef¬ 
ficiency  and  attenuation  were  checked  against  hand  calcula¬ 
tions.  Scatter  was  not  considered. 

Projection  images  were  simulated  for  step- and- shoot 
mode  every  3°  over  a  half-circular  trajectory  by  rotating  the 
detector  from  right  lateral  over  the  chest  to  left  lateral  as 
shown  in  Fig.  4.  Bins  in  projection  images  were  0.36  cm 
wide.  Finite  bin  width  was  modeled  by  tracing  four  cones 
from  each  bin,  i.e.,  one  cone  from  the  center  of  the  four 
quadrants;  thus  in  total  4X  121=484  rays  were  traced  from 


Fig.  4.  Detector  trajectory  spanning  180°  anteriorly  is  displayed  with  over¬ 
lapping  gray  bars.  Axis  of  rotation  is  marked  with  an  asterisk.  Common 
volume — voxels  projected  onto  the  detector  at  every  view — is  inside  the 
dashed  line. 


each  detector  bin.  The  detector-to-image-plane  radius  of  ro¬ 
tation  was  26.5  cm  with  the  axis  of  rotation  centered  near  the 
phantom  spine.  Using  this  detector  and  trajectory,  there  was 
no  truncation  of  the  transaxial  slices  containing  tumors,  as 
demonstrated  in  Fig.  4  by  the  dashed  line  that  encloses  the 
common  volume,  i.e.,  the  region  that  is  projected  onto  the 
detector  surface  at  every  angle.  These  simulated  noise-free 
2D  projection  images  were  degraded  with  varying  levels  of 
Poisson  noise  corresponding  to  scan  times  of  4,  8,  and  20 
min.  Ensembles  of  25  independent  noisy  projection  images 
were  generated  for  each  combination  of  scan  time  and  tumor 
size. 

Image  ensembles  were  reconstructed  via  ordered- subsets 
expectation  maximization  (OSEM)  (Ref.  27)  using  six  sub¬ 
sets  and  five  iterations,  which  are  in  the  range  of  typically 
used  values.  Attenuation  correction  was  achieved  by  mod¬ 
eling  photon  detection  probabilities  within  OSEM.  In  the  re¬ 
constructed  images,  the  intensities  of  deep  background  re¬ 
gions  were  similar  to  those  in  the  phantom,  thus  validating 
attenuation  correction.  Our  reconstructions  did  not  model 
distance-dependent  spatial  resolution,  as  is  typical  of  most 
clinical  SPECT  imaging.  Image  voxels  were  0.36  cm  wide,  a 
factor  of  2  coarser  than  voxels  of  grids  on  which  phantoms 
were  placed  for  simulating  projections.  During  reconstruc¬ 
tion,  four  rays  were  traced  from  every  projection  bin  through 
the  reconstruction  grid  in  order  to  average  over  the  finite 
detector  bin  width.  Rays  were  oriented  perpendicular  to  the 
image  plane  and  spaced  equally  from  one  another  in  the  tran¬ 
saxial  dimension  of  a  bin.  Additional  rays  were  not  traced 
along  the  axial  (superior-inferior)  dimension  because  of  the 
exact  axial  alignment  between  projection  bins  and  image 
voxels  at  every  detector  view.  Noisy  reconstructed  images 
were  smoothed  with  3D  Gaussian  kernels  of  8,  14,  and  20 
mm  FWHMs. 


II.C.  Image  analysis 

Contrast-to-noise  ratios  (CNRs)  and  localization  metrics 
were  calculated  from  noisy  image  ensembles.  For  CNR  cal¬ 
culations,  3D  ROIs  were  drawn  around  the  tumor.  Approxi¬ 
mately  7.2  mm  away  from  these  ROI  boundaries,  shell-like 
ROIs  were  defined  in  the  background.  ROIs  were  approxi¬ 
mately  spherical  for  stationary  tumors.  For  lung  tumors, 
ROIs  were  irregularly  shaped  so  as  to  account  for  the  asym¬ 
metric  blurring  effects  of  respiratory  motion.  Both  tumor  and 
background  ROIs  were  the  size  of  tumor  volumes  and  en¬ 
compassed  21,  45,  and  139  voxels  for  respective  tumor  di¬ 
ameters  of  10.8,  14.4,  and  21.6  mm. 

Contrast  was  defined  as 


where  Tt  and  Bi  are  mean  tumor  and  background  ROI  inten¬ 
sities  in  image  i.  Mean  contrast, 
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C  = 


y\N  r 

Zji=V^i 


N 


(2) 


was  calculated  from  the  N=  25  noisy  reconstructed  images. 
This  definition  of  contrast  allows  for  negative  values  when 
the  background  is  hotter  than  the  tumor  in  order  to  avoid 
upward  bias  of  mean  contrast  values  when  background  noise 
dominates  the  tumor  signal.  Noise  was  calculated  as  the  stan¬ 
dard  deviation  of  contrast: 


(Tc  = 


s£t(  Q-C)2 

N-  1 


(3) 


Ensemble  average  contrast-to-noise  ratio  was 
C 

CNR  =  — . 

°c 


(4) 


Localization  accuracy  and  precision  were  assessed  by  cal¬ 
culating  the  mean  and  standard  deviation  of  distances  be¬ 
tween  true  and  measured  tumor  centroids  over  ensembles  of 
noisy  reconstructed  images.  Tumor  locations  in  noisy  images 
were  selected  via  a  forced  choice  task  using  nonprewhitening 
filters.  This  approach  is  relevant  to  radiation  therapy  de¬ 
livery  because  disease  is  known  to  exist;  yet  the  exact  posi¬ 
tion  varies  with  each  time  the  patient  is  positioned  on  the 
treatment  table. 

Since  the  approximate  tumor  location  is  known  from 
treatment-planning  images,  the  search  can  be  limited  to  a 
relatively  small  volume.  In  this  study,  3D  search  volumes — 
14.4-mm-radius  spherical  ROIs — were  defined  around  each 
tumor.  The  size  of  each  search  volume  was  intended  to  en¬ 
compass  the  uncertainties  in  tumor  location  that  are  typical 
for  radiation  therapy.  Each  search  volume  contained  257 
voxels  that  were  subdivided  three  times  in  each  dimension. 
The  location  in  the  ROI  that  had  the  highest  cross  correlation 
(XC)  with  a  noise-free  tumor  template  was  selected  as  the 
measured  tumor  centroid: 


XC,  = 


template.  •  image/ 


V(template.-  •  template.) (image/  •  image/) 


(5) 


r m,i =  arg  max(XC/j) ,  (6) 

j^n 

where  template^  represents  the  template  centered  at  voxel  j 
within  ROI  12  of  ensemble  image  i.  The  denominator  in  the 
above  cross-correlation  equation  normalizes  the  response 
such  that  it  does  not  depend  on  voxel  intensities.  As  such 
cross-correlation  values  were  suppressed  for  hot  background 
structures  (e.g.,  heart)  that  would  have  otherwise  elicited 
greater  responses  than  the  nearby  tumor.  Tumor  templates 
were  generated  by  convolving  true  tumor  geometries  with  a 
14.4-mm-Gaussian  kernel  that  approximates  typical  spatial 
resolutions  in  these  images.  For  lung  tumors,  the  true  tumor 
geometry  included  asymmetric  blurring  from  respiratory  mo¬ 
tion.  Tumor  templates  were  constructed  on  cubic  grids  of 
0.12-cm-wide  voxels  with  grid  widths  of  27,  33,  and  39  vox¬ 
els  for  tumor  diameters  of  10.8,  14.4,  and  21.6  mm,  respec¬ 
tively. 
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Fig.  5.  Sample  reconstructed  images,  smoothed  with  14-mm-FWHM 
Gaussian,  for  each  combination  of  tumor  size  and  scan  time. 


Mean  localization  error  dr  was  defined  as  the  offset  dis¬ 
tance  dr  between  true  rt=(xt,yt,zt)  and  measured  tumor  cen¬ 
troids  rmi=(xmi,ymi,zmj)  averaged  over  the  ensemble  of  N 
=  25  noisy  reconstructed  images: 

dr  =  V(xm  -  x()2  +  (ym  -  y ()2  +  (zm  -  z,)2,  (7) 


_  l 

dr=  ~2j  dr:. 


(8) 


Standard  deviation  in  localization  error  was  computed  as 


® dr 


2^  (Jr/-  dr)2 

N-  1 


(9) 


III.  RESULTS 

III.A.  Reconstructed  images 

Reconstructed  images  for  each  combination  of  tumor  di¬ 
ameter  and  scan  time  are  presented  in  Fig.  5.  These  images 
are  noisy  realizations,  not  ensemble  averages,  and  have  been 
smoothed  with  a  3D,  14-mm-FWHM  Gaussian  kernel.  Grids 
are  superimposed  on  images  such  that  each  square  element  is 
centered  on  a  tumor. 

Qualitative  trends  are  apparent  in  these  images.  As  scan 
time  increases,  noise  decreases,  thereby  improving  the  tumor 
visibility.  For  instance,  with  the  4  min  scans,  image  noise 
resembles  the  expected  signal  for  10.8  and  14.4  mm  diameter 
tumors  in  some  locations,  especially  in  the  lungs  and  near 
posterior  of  the  phantom.  When  scan  time  is  increased  to  8 
and  then  20  min,  there  are  fewer  and  fewer  hot  spots  off 
center  of  the  grid,  which  suggests  that  the  apparent  activity  is 
the  true  tumor  signal.  Tumor-to-background  uptake  ratios 
have  a  substantial  effect  on  tumor  visibility,  as  those  on  the 
patient  right  side  with  6:1  uptake  ratios  are  more  noticeable 
than  left-sided  tumors  with  3:1  uptake  ratios.  Also  note  the 
impact  of  tumor  location  relative  to  the  detector  trajectory: 
Tumors  near  the  chest  wall,  which  are  on  average  closer  to 
the  detector,  are  strikingly  more  apparent  than  lung  tumors 
and  stationary  tumors  near  the  spine.  The  detector  trajectory 
also  affects  apparent  tumor  shape.  Anterior  tumors  are  ap- 
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Fig.  6.  CNRs  calculated  from  en¬ 
sembles  of  noisy  reconstructed  images 
and  presented  as  a  function  of  location 
for  each  combination  of  scan  time  and 
tumor  diameter. 


proximately  circular,  while  other  visible  tumors  are  elliptical 
with  varying  degrees  of  eccentricity.  Greater  eccentricity  oc¬ 
curs  when  the  detector-to-tumor  distance  is  substantially  dif¬ 
ferent  over  the  detector  trajectory. 

III.B.  Contrast-to-noise  ratio 

CNRs  were  calculated  from  ensembles  of  noisy  recon¬ 
structed  images  using  Eqs.  (l)-(4).  CNR  values  are  dis¬ 
played  in  Fig.  6  as  a  function  of  anatomy  for  the  investigated 
scan  times  and  tumor  diameters.  Many  of  the  visual  trends 
noted  in  the  reconstructed  images  from  Fig.  5  are  expressed 
quantitatively  in  Fig.  6  as  ensemble  CNRs.  CNRs  increase 
approximately  by  the  square  root  of  imaging  time.  Right¬ 
sided  tumors  with  6:1  radiotracer  uptake  ratios  have  higher 
CNRs  than  their  laterally  symmetric  counterparts  of  3:1  up¬ 
take  by  an  average  factor  of  1.6.  CNRs  of  anterior  tumors  are 
typically  3.5  times  greater  than  their  posterior  counterparts. 
These  trends  are  shown  in  Fig.  7,  thus  demonstrating  the 


0  2  4  6 

cnr 

Fig.  7.  Differences  in  CNR  due  to  radiotracer  uptake  and  anatomical  loca¬ 
tion.  Circles  represent  differences  in  CNR  between  anterior  tumors  and  their 
posterior  counterparts.  Diamonds  indicate  the  CNRs  of  laterally  symmetric 
tumors  that  differ  by  radiotracer  uptake. 


relative  impact  of  radiotracer  uptake  versus  the  detector  tra¬ 
jectory  on  CNR. 

III.C.  Localization 

Focalization  metrics — the  mean  localization  error  dr  and 
ensemble  standard  deviation  crdr — were  calculated  using  Eqs. 
(5)— (9)  and  are  reported  in  Fig.  8.  Mean  errors  are  repre¬ 
sented  by  white  and  black  boxes  that  respectively  correspond 
to  tumors  of  6:1  and  3:1  radiotracer  uptakes.  Error  bars  in¬ 
dicate  standard  deviations.  The  different  tumor  sizes  are 
separated  by  rows.  Columns  differentiate  tumor  locations. 
Note  the  diagram  above  each  column  in  Fig.  8  that  specifies 
anatomical  locations.  Within  each  column  there  are  six  data 
points — three  white  and  three  black  boxes — that  correspond 
to  scan  times  of  4,  8,  and  20  min  for  each  uptake  ratio.  When 
viewing  Fig.  8  from  left  to  right,  tumor  locations  vary  from 
anterior  to  lateral  to  posterior.  A  bold  dashed  line  separates 
stationary  and  lung  tumors.  The  relative  order  of  lung  tumors 
is  the  same  as  that  of  stationary  tumors. 

In  all  but  the  regions  of  poorest  SPECT  image  quality, 
tumors  of  6:1  uptake  are  typically  localized  as  accurately 
with  4  min  of  scan  time  as  tumors  of  3:1  uptake  that  are 
imaged  for  20  min.  This  trend  is  demonstrated  in  Fig.  8  by 
observing  the  relative  height  of  boxes,  specifically  the  first 
white  and  the  third  black  boxes  within  a  column.  For  loca¬ 
tions  at  which  noise  and  blur  are  severe,  localization  is  de¬ 
termined  primarily  by  random  chance  for  both  3:1  and  6:1 
uptake  ratios.  In  these  cases,  localization  accuracy  is  similar 
and  poor  for  the  two  uptake  ratios.  This  effect  is  illustrated 
by  the  10.8  mm  tumors  in  the  lungs. 

Tumor  location  has  a  major  impact  on  localization  accu¬ 
racy.  In  Fig.  8,  there  is  a  general  trend  where  localization 
errors  worsen  when  tumor  location  changes  from  anterior  to 
lateral  to  posterior.  As  with  CNRs,  tumor  location  relative  to 
the  detector  trajectory  can  have  a  larger  effect  on  localization 
error  than  radiotracer  uptake  ratio.  These  effects  are  shown 
in  Fig.  8  when  making  comparisons  among  stationary  tumors 
near  the  chest  wall  of  different  uptake  ratios,  and  then  with 
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Fig.  8.  Ensemble  localization  errors  for  different  combinations  of  tumor  size,  location,  radiotracer  uptake,  and  imaging  time.  Mean  errors  are  displayed  as 
black  and  white  boxes  for  respective  uptake  ratios  of  3:1  and  6:1.  Error  bars  indicate  standard  deviations.  Tumor  sizes  are  separated  by  rows  and  tumor 
locations  by  columns.  Within  each  column,  there  are  six  data  points  that  correspond  with  scan  times  of  4,  8,  and  20  min  for  laterally  symmetric  tumors  of 
different  uptake  ratios. 


counterparts  near  the  spine  and  posterior  ribs.  Though  local¬ 
ization  is  generally  poor  for  posterior  tumors,  certain  results 
are  encouraging.  When  using  relatively  short  4  min  scans, 
anterior  tumors  with  diameters  of  14.4  and  21.6  mm  and  6:1 
uptake  ratios  had  mean  localization  errors  less  than  2  mm. 

III.D.  Localization  accuracy  as  a  function  of 
contrast-to-noise  ratio 

Mean  localization  errors  are  plotted  as  a  function  of 
CNRs  in  Fig.  9.  There  is  a  noticeable  trend.  When  CNRs  are 
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Fig.  9.  Mean  localization  errors  plotted  as  a  function  of  CNRs.  Tumors  of 
3:1  and  6:1  uptakes  are  differentiated  by  square  and  diamond  markers. 
Markers  that  are  crossed  signify  lung  tumors,  while  solid  markers  represent 
stationary  tumors. 


close  to  zero,  localization  errors  are  around  11  mm,  which  is 
the  expected  error  when  randomly  selecting  radial  locations 
from  a  14.4-mm-radius  search  volume.  In  this  regime,  noise 
dominates  and  SPECT  provides  very  little  information.  As 
CNRs  increase  to  around  5,  localization  errors  improve  rap¬ 
idly  because  SPECT  is  providing  more  information.  Im¬ 
provements  in  localization  are  more  modest  as  CNRs  con¬ 
tinue  to  increase.  In  this  region,  there  is  comparatively  less  to 
gain  from  enhancing  the  already  pronounced  tumor  signal. 

There  are  other  more  subtle  effects  in  Fig.  9.  In  the  CNR 
range  of  3-5,  localization  is  typically  better  for  tumors  of  6:1 
than  3:1  uptake.  Localization  error  is  generally  better  for 
stationary  tumor  of  the  same  uptake  ratio  than  lung  tumors 
with  similar  CNRs.  Two  stationary  tumors  of  3:1  uptake  with 
relatively  high  CNRs  are  noticeably  skewed  above  the  other 
data  points.  The  unexpected  degradation  in  localization  error 
is  attributable  to  proximity  to  the  heart  where  radiotracer 
concentrations  are  much  higher  and  suppress  tumor  signal 
inferiorly.  These  results  show  potential  challenges  for  local¬ 
izing  targets  near  a  hot  background. 

IV.  DISCUSSION 

In  this  simulation  study,  on-board  SPECT  was  investi¬ 
gated  for  localizing  functional  and  molecular  targets.  Several 
parameters  were  considered:  Scan  time  and  tumor  location, 
diameter,  and  uptake  ratio.  Among  these,  one  of  the  most 
influential  factors  was  radiotracer  uptake  in  tumor  relative  to 
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background:  As  indicated  by  CNRs,  localization  metrics,  and 
visual  inspection.  Tumors  of  6:1  uptake  were  strikingly  more 
noticeable  than  their  laterally  symmetric  counterparts  of  3:1 
uptake.  These  visual  observations  were  supported  by  quanti¬ 
tative  values — CNRs  and  localization  accuracy.  CNRs  were 
on  average  1.6  times  greater  for  the  tumors  of  6:1  radiotracer 
specificity.  The  observed  ratio  of  1.6  is  less  than  the  true 
phantom  contrast  ratio  of  2.5.  In  the  absence  of  spatial  reso¬ 
lution  modeling,  blurring  from  limited  spatial  resolution  has 
a  relatively  greater  impact  on  contrast  recovery  of  6:1  tumors 
compared  with  those  of  3:1  uptake.  This  result  is  consistent 
with  the  findings  of  Ref.  32  where  scatter,  a  type  of  blurring, 
more  substantially  reduced  the  recovered  contrast  of  hotter 
lesions. 

Typically,  tumors  of  6:1  uptake  were  localized  as  accu¬ 
rately  using  4  min  scans  as  3:1  tumors  that  had  been  imaged 
for  20  min.  This  difference  in  time  has  important  implica¬ 
tions,  especially  in  the  context  of  radiation  therapy  delivery 
where  short  scans  are  important  in  minimizing  patient  mo¬ 
tion  and  ensuring  machine  throughput.  Thus,  radiotracer 
specificity  will  be  an  important  criterion  when  selecting  or 
developing  potential  radiotracers  for  on-board  SPECT  imag¬ 
ing.  Conversely,  diagnostic  or  treatment-planning  SPECT 
images  could  be  used  to  estimate  tumor  characteristics — 
radiotracer  uptake  and  size — that  affect  localization  in  a  pre¬ 
dictable  way  in  order  to  gauge  whether  individual  patients 
would  benefit  from  on-board  SPECT.  Further,  diagnostic  or 
treatment-planning  SPECT  images  could  provide  radiotracer- 
concentration  templates  for  estimating  localization.  Such 
templates  may  be  particularly  useful  for  localizing  tumors 
with  inhomogeneous  radiotracer  uptake  that  is  comparable  in 
width  to,  or  wider  than,  the  on-board  SPECT  spatial  reso¬ 
lution. 

This  study  utilized  a  180°  orbit  in  order  to  avoid  viewing 
the  patient  through  the  couch.  Radiation  therapy  couches  are 
thicker  and  wider  than  typical  diagnostic  imaging  couches, 
such  that  SPECT  data  acquired  through  the  couch  would  be 
substantially  attenuated,  and  the  SPECT  detector  would  be 
displaced  away  from  the  patient,  thereby  degrading  spatial 
resolution.  Averaging  over  all  angles  of  the  180°  orbit,  tu¬ 
mors  near  the  chest  wall  are  on  average  closer  to  the  detec¬ 
tor,  which  implies  better  spatial  resolution,  and  they  are  less 
attenuated  by  tissue.  Because  of  attenuation,  image  noise 
was  greater  near  the  spine  than  the  chest  wall.  Quantum 
noise  is  not  reversible  even  with  attenuation  correction.  As  is 
typical  for  most  clinical  SPECT  imaging,  the  reconstructions 
did  not  model  spatial  resolution.  Consequently,  poorer  spatial 
resolution  and  poorer  contrast  are  expected  for  the  posterior 
tumors.  Modeling  spatial  resolution  would  increase  image 
reconstruction  time — an  important  consideration  for  on¬ 
board  SPECT.  That  noted,  spatial  resolution  modeling  is 
worth  evaluating  in  a  future  study,  since  it  has  potential  to 
improve  localization.3'  Because  localization  is  better  in  re¬ 
gions  proximal  to  the  detector  trajectory  and  since  tumor 
location  is  known  approximately,  on-board  SPECT  trajecto¬ 
ries  could  be  optimized  for  imaging  specific  tumor  sites. 

CNR  and  localization  error  are  related,  but  one  metric 
may  not  necessarily  be  a  good  predictor  of  the  other.  In  this 


study,  noisy  reconstructed  images  were  smoothed  using 
Gaussian  kernels  with  FWHMs  of  8,  14,  and  20  mm.  CNRs 
and  localization  errors  were  calculated  for  noisy  images  and 
for  images  with  each  degree  of  smoothing.  CNRs  improved 
with  smoothing.  In  contrast,  localization  errors  improved  at 
certain  anatomical  locations  but  worsened  for  others.  On  av¬ 
erage  localization  errors  were  comparable  across  different 
degrees  of  smoothing.  These  results  demonstrate  that  local¬ 
ization  errors  cannot  be  estimated  directly  from  CNRs.  En¬ 
couragingly,  these  results  also  suggest  that  localization  accu¬ 
racy  is  somewhat  robust  to  smoothing.  It  follows  that  they 
may  be  robust  to  iteration  number  when  using  iterative  re¬ 
construction  since  smoothness  is  highly  correlated  with  itera¬ 
tion  number.  Smoothing  did,  however,  impact  localization 
for  certain  tumors,  and  these  effects  need  to  be  studied  fur¬ 
ther. 

CT  provides  anatomical  information  that  would  be  valu¬ 
able  for  correcting  nonuniform  attenuation  in  on-board 
SPECT  images.  Attenuation  coefficients  could  be  estimated 
from  on-board  cone  beam  CT  or  from  registered  CT  images. 

We  do  not  anticipate  that  on-board  SPECT  imaging  would 
be  used  for  every  fraction  in  a  highly  fractionated  treatment 
because  of  associated  imaging  time  and  patient  dose.  On¬ 
board  SPECT  is  better  suited  for  treatments  with  a  limited 
number  of  fractions  such  as  stereotactic  body  radiation 
therapy  (SBRT).  SBRT  is  used  to  treat  a  variety  of  tumors, 
including  breast  oligometastasis. 34,35  For  SBRT,  imaging 
time  is  less  of  an  issue  because  additional  time  is  allotted  for 
patient  positioning  as  compared  with  most  highly  fractioned 
treatments.  Concerning  patient  dose,  a  typical  25  mCi  injec¬ 
tion  of  "mTc-sestamibi36  results  in  an  effective  dose  of  7.9 
mSv,37  which  is  about  1/3  the  effective  dose  of  an  on-board 
CBCT  chest  scan.  Cost  is  also  an  issue  that  would  have  to 
be  evaluated  along  with  any  improvements  in  morbidity 
and/or  mortality  through  clinical  studies. 

Localization  values  may  change  depending  on  the  coarse¬ 
ness  of  image  voxels  or  the  size  of  the  search  volume.  In  this 
study,  centroids  were  estimated  on  a  voxelized  grid,  limited 
by  the  half  width  of  a  sub  voxel.  Localization  would  likely 
worsen  by  using  a  larger  search  volume  with  more  locations 
where  random  noise  fluctuations  could  mimic  tumor,  particu¬ 
larly  for  barely  visible  tumors. 

Our  simulation  study  makes  absolute  estimates  regarding 
localization.  Localization  values  will  change — likely  for  the 
worse  as  additional  factors  such  as  scatter  (estimated  to  ac¬ 
count  for  20%-40%  of  the  detected  counts  in  scanner 
data39),  tumor  diffuseness,  and  tumor  inhomogeneities  are 
considered.  Localization  accuracy,  however,  will  likely 
change  for  the  better  as  trajectories  are  focused  on  the  target 
region,  as  scatter  correction  is  implemented,  and  as  estima¬ 
tion  methods  are  developed  specifically  for  tumor  inhomoge¬ 
neities.  This  current  study  contributes  as  a  significant  first 
step,  by  incorporating  a  significant  degree  of  realism  (e.g., 
spatial  resolution  and  attenuation),  clearly  defining  that  level 
of  realism,  and  then  returning  absolute  numbers  for  localiza¬ 
tion  accuracy  and  precision.  Such  absolute  estimates  are  nec¬ 
essary  in  order  to  access  whether  on-board  SPECT  is  even 
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possibly  feasible,  and  these  estimates  form  a  baseline  which 
will  help  in  quantifying  and  understanding  the  effects  of  fu¬ 
ture  developments. 

V.  CONCLUSION 

SPECT  imaging  was  investigated  for  localizing  functional 
and  molecular  targets  immediately  prior  to  radiation  therapy. 
CNRs  and  localization  errors  were  analyzed  as  a  function  of 
scan  time  and  tumor  size,  location,  and  radiotracer  uptake 
using  computer  simulations.  Localization  errors  were  less 
than  2  mm  for  certain  tumors  using  relatively  short  4  min 
scans.  These  encouraging  results  warrant  further  investiga¬ 
tion  of  on-board  SPECT  for  localizing  functional  and  mo¬ 
lecular  targets. 
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ABSTRACT 


Purpose:  Single  photon  emission  computed  tomography  (SPECT)  is  being  investigated 
for  imaging  inside  radiation  therapy  treatment  rooms  to  localize  biological  targets.  Here 
computer  simulations  were  used  to  analyze  locational  and  directional  dependencies  in 
localization  errors  and  to  assess  the  effects  of  spatial  resolution  modeling  and  observer 
40  normalization  on  localization  performance. 

Methods:  SPECT  images  of  the  XCAT  phantom,  containing  12  hot  tumors,  were 
reconstructed  with  detector  response  function  compensation  (DRC)  and  without  DRC 
(nDRC).  Numerical  observers  were  forced  to  select  the  most  suspicious  tumor  location  - 
using  normalized  (NXC)  or  un-normalized  (XC)  cross  correlation  -  from  3-cm-diameter 
45  search  volumes  that  each  contained  only  one  tumor.  For  each  tumor  site,  localization  was 
optimized  as  a  function  of  iteration  number  and  post-reconstruction  smoothing. 
Localization  error  -  the  distance  between  true  and  estimated  tumor  positions  -  was 
calculated  across  the  ensembles  of  80  images.  Direction-dependent  localization  bias  and 
precision  were  estimated  from  the  image  ensemble. 

50  Results:  For  the  6  superficial  tumors  in  close  proximity  to  the  detector  trajectory,  mean 
localization  errors  were  <  2  mm  and  were  lowest  or  comparable  using  DRC -NXC,  though 
differences  from  DRC-XC  and  nDRC -NXC  were  not  statistically  significant.  DRC -NXC 
did  provide  statistically  significantly  better  localization  than  nDRC-XC  for  5  of  these  6 
tumors.  At  the  other  6  sites  where  attenuation  was  more  severe  and  the  distance 
55  generally  greater  between  tumor  and  detector,  DRC  typically  did  not  show  better 

localization  than  nDRC.  Observer  normalization  improved  localization  substantially  for 
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a  tumor  near  the  hotter  heart.  Localization  errors  were  anisotropic  and  dependent  on 
tumor  location  relative  to  the  detector  trajectory. 

Conclusions:  This  computer-simulation  study  compared  localization  performance  for 
60  normalized  and  un-normalized  numerical  observers,  which  were  used  to  estimate  tumor 
positions  in  SPECT  images,  reconstructed  with  and  without  DRC.  For  tumors  localized 
to  <  2  mm  on  average,  which  are  good  candidates  for  SPECT-guided  radiation  therapy, 
localization  performance  typically  improved  by  compensating  for  the  detector  response 
function  and  by  using  a  normalized  observer.  The  observed  direction-dependent 
65  localization  errors  have  important  implications  for  radiation  therapy  and  are  relevant  to 
SPECT  imaging  in  general. 

Keywords:  SPECT,  localization,  direction-dependent  localization  errors,  detector 
response  function  compensation,  numerical  observer 
I.  INTRODUCTION 

70  Biologically-conformal  radiation  therapy  has  been  proposed  for  sculpting  dose  in 

regions  of  high  tumor  burden,  increased  growth,  and  hypoxia.1  Currently,  however,  there 
are  no  imaging  modalities  inside  radiation  therapy  treatment  rooms  that  provide 
functional  information  for  localizing  biological  targets.  In  a  previous  study,  we 
investigated  certain  aspects  of  the  technical  feasibility  of  using  single  photon  emission 
75  computed  tomography  (SPECT)  for  the  task  of  localizing  biological  targets  while  a 
patient  is  on  the  treatment  couch,  immediately  prior  to  radiation  delivery.2  In  this 
previous  study,  SPECT  images  were  reconstructed  without  compensation  for  the  detector 
response  function,  and  normalized  numerical  observers  were  used  to  assess  localization 
errors  as  a  function  of  the  location,  size,  and  uptake  ratio  of  hot  tumors. 


80  In  the  present  study,  the  evaluation  of  localization  errors  was  extended  to 

directional  dependencies,  which  were  analyzed  relative  to  the  detector  trajectory  and  the 
attenuating  medium.  Also  evaluated  for  localization  performance  were  the  use  of  spatial 
resolution  modeling  within  image  reconstruction  and  the  use  of  normalization  in  the 
numerical  observer.  Assessments  of  direction-dependent  localization  errors  are  important 
85  to  external  beam  radiation  therapy  because  radiation  is  typically  delivered  from  multiple 
gantry  angles. 

Related  work  on  spatial  resolution  modeling  includes  localization  receiver 
operating  characteristic  (LROC)  studies.3  These  LROC  studies  involved  human 

rn  o 

observers  and  computer  simulations  of  Ga  SPECT  tumor  imaging  using  a  360  detector 
90  trajectory.  Observers  considered  a  broad  region  of  the  mediastinum  where  lymphoma 
commonly  occurs,  and  LROC  performance  was  averaged  over  many  tumor  sites.  These 
studies  found  improved  LROC  performance  when  utilizing  SPECT  images  reconstructed 
with  detector  response  compensation  (DRC),  as  compared  to  utilizing  images  with  no 
DRC  (nDRC). 

95  II.  METHODS  AND  MATERIALS 

A.  Phantom 

A  female  XCAT  Phantom,  an  extended  version  of  the  NCAT  Phantom4,  was 
implemented  on  a  grid  of  0.15  x  0.15  x  0.15  cm3  voxels.  Twelve,  1.5-cm-diameter 
tumors  were  simulated  -  in  axilla,  breast,  chest  wall,  lung,  rib,  vertebra,  and  other  soft 
100  tissues  -  for  a  broad  sampling  of  locations  deep,  superficial,  anterior,  and  posterior  in  an 
axial  slice  superior  to  the  heart.  Tumor  activity  concentration  was  1.5  pCi/cnr3,  based  on 
breast  tumor  imaging  with  "mTc  cardiac  perfusion  radiotracers.5  Relative  activity  ratios 


were  20:6:1  for  heart,  tumor,  and  other  tissues.5, 6  The  simulated  activity  distribution  and 
the  12  tumor  sites  labeled  A-L  are  shown  in  Figure  1.  Beneath  the  supine  phantom  with 
105  raised  arms,  a  flat-top  couch  was  simulated  that  had  a  lateral  width  of  50.5  cm. 

Attenuation  was  modeled  for  140  keV  photons,  the  emission  energy  of  "mTc.  Linear 
attenuation  coefficients  were  as  follows  for  the  different  structures:  lung,  0.043  cm'1;  soft 
tissue,  0.15  cm'1;  rib,  0.22  cm'1;  vertebra,  0.18  cm'1,  and  couch  0.25  cm'1.  Respiratory 
and  cardiac  motions  were  not  modeled. 

110  B.  Image  Simulation 

In  SPECT,  a  collimator  is  used  to  select  photons  emitted  in  certain  directions. 

The  low-energy,  high-resolution  (LEHR)  collimator  simulated  in  this  study  had  circular 
parallel  holes,  0.14  cm  in  diameter  and  2.635  cm  in  effective  length.  Hole  length  /  and 
diameter  d  determine  angles  -  <p(-tanA(d/[),  tan'1  (d/l))  -  through  which  rays  may  pass 
115  directly  to  the  detector.  Rays  perpendicular  to  the  detector  surface,  where  (f)  is  zero,  have 
the  highest  probability  of  being  accepted.  Angular  ray  weights  a((/))  for  circular  parallel 
holes  are  described  by 


2  cos  1 


/•tan  (/> 


/•tan  (f) 


1- 


/•tan  (f) 


(1) 


where  c  is  a  normalization  factor.7 

Figure  2  shows  schematics  of  (a)  the  collimator  hole  geometry  and  (b)  the  cone  of 
120  rays  used  to  simulate  distance-dependent  spatial  resolution.  The  cone  of  rays  was  divided 
into  49  segments  of  equal  area.  Rays  were  weighted  by  the  mean  a((f))  of  each  segment. 


Intrinsic  blur  was  modeled  with  a  2D,  0.32-cm-FWHM  Gaussian  kernel.  Figure  2c 
shows  the  agreement  between  measured  and  simulated  line  source  profiles  at  distances  of 


15  to  30  cm.  Profiles  -  1 ,8-cm  thick  -  were  drawn  perpendicular  to  the  long  axis  of 
125  simulated  and  physical  line  sources  that  were  4  cm  in  length  and  0.1  cm  in  diameter.  The 
physical  source,  a  capillary  tube,  was  filled  with  "mTc.  Measurements  were  made  on  a 
Trionix  Triad  scanner,  equipped  with  parallel-hole  LEHR  collimation.  Projection  images 
were  generated,  with  detector  resolution  modeled  as  just  described,  using  the  software 
SPECT-MAP,  which  analytically  models  collimator  efficiency,  detector  efficiency, 

130  radioactive  decay,  isotope  branching,  non-uniform  attenuation,  and  distance-dependent 
spatial  resolution.  Scatter  was  not  modeled  in  the  phantom  or  detector.  Septal 
penetration  was  not  simulated.  SPECT-MAP  has  been  widely  utilized  and  tested.2, 8-12 

A  conformal  detector  trajectory  was  defined  such  that  for  each  detector  view,  the 
detector  radius  was  adjusted  to  achieve  0.5-cm  proximity  between  the  patient-side  surface 
135  of  the  collimator  and  the  phantom  or  couch.  Detector  axial  width  was  considered  when 
varying  the  detector  radius,  thereby  avoiding  collisions  with  breast,  raised  arms,  and 
couch.  As  measured  from  the  detector  axis  of  rotation  near  the  spine,  the  distance  to  the 
collimator  surface  ranged  from  a  minimum  of  23.9  cm  in  the  anterior  above  the  breasts  to 
a  maximum  of  27.4  cm  near  the  arms.  There  were  60  detector  views  over  the  anterior 
140  180°  trajectory,  as  shown  in  Figure  3.  The  simulated  54-cm  detector  width  was 

sufficiently  broad  as  to  avoid  image  truncation  transaxially.  In  the  axial  dimension,  the 
simulated  detector  was  more  compact  at  15  cm  than  a  full-size  gamma  camera  of  40  cm 
width.13  Projection  bins  were  0.15  x  0.15  cm2.  Noise-free  projection  images  were 
degraded  with  Poisson  noise  to  simulate  5  seconds  of  imaging  time  per  detector  view.  In 
145  regions  of  the  detector  illuminated  by  background,  i.e.,  not  by  heart  or  tumor, 


approximately  8  counts  per  second  were  recorded  per  cm2.  An  ensemble  of  80  noisy 
realizations  was  generated  using  this  detector  trajectory. 

C.  Image  Reconstruction 

Images  were  reconstructed  using  OSEM  14  and  then  smoothed  with  a  3D  Gaussian 
150  filter.  For  each  tumor  location  and  each  combination  of  DRC  versus  nDRC  and 

normalized  versus  un-normalized  observer,  localization  was  evaluated  across  a  broad 
range  of  image  updates  (i.e.  products  of  numbers  of  iterations  and  subsets)  and  of 
smoothing  kernel  widths,  in  order  to  determine  an  optimal  number  of  updates  and  degree 
of  smoothing  at  each  location  and  for  each  method  of  reconstruction  and  each  method  of 
155  localization  estimation.  These  different  methods  were  compared  using  the  iteration 

numbers  and  smoothing  widths  that  were  optimal  for  each  method.  In  the  comparison 
results  presented  below,  four  subsets  were  used  to  accelerate  DRC  reconstructions,  but 
only  one  subset  was  used  for  nDRC  images  because  the  iteration-number  optimization 
indicated  that  1  iteration  with  4  subsets  gave  too  many  image  updates  for  optimal 
160  localization  at  several  tumor  sites.  Note  that  one-subset  OSEM  is  MLEM.  1416  Post¬ 
reconstruction  smoothing  was  performed  using  3D  Gaussian  kernels  that  ranged  in 
FWHM  from  0  to  2.5  cm  in  increments  of  0.5  cm.  In  total,  80  noisy  images  were 
reconstructed  using  2  variations  of  OSEM  (DRC  and  nDRC)  and  smoothed  to  different 
degrees  with  6  Gaussian  kernels  at  iterations  1-25. 

165  Both  reconstruction  methods  used  the  same  ensemble  of  noisy  projection  images. 

As  such,  reconstructed  images  were  paired  with  the  only  difference  being  the 
reconstruction  method  (DRC  or  nDRC).  For  nDRC,  central  rays  were  traced 
perpendicular  to  the  image  plane  and  from  each  detector  bin  through  the  phantom.  The 


170 


DRC  model  accounted  for  distance-dependent  collimator  spatial  resolution  and  crystal 
blurring  and  was  matched  exactly  to  that  used  for  simulating  projection  images.  For  both 
reconstruction  methods,  non-uniform  attenuation  was  calculated  along  each  ray.  The 
attenuation  map  was  ideal.  Reconstructed  images  were  implemented  on  a  grid  of  0.30  x 
0.30  x  0.30  cm3  voxels. 

Tumor  positions  were  estimated  in  background-subtracted  noise-free  images,  as  a 
175  method  for  verifying  that  tumor  positions  in  reconstructed  images  were  consistent  with 
those  in  the  phantom.  Noise-free  nDRC  and  DRC  images  were  reconstructed  using  the 
methods  described  above.  Noise-free  images  without  tumor  present  were  used  to 
estimate  background  activity.  In  background-subtracted  noise-free  images,  1.8-cm- 
diameter  ROIs  were  placed  on  the  1 ,5-cm-diameter  tumors.  The  center  of  mass  was 
180  calculated  for  each  tumor  site  and  across  the  25  iterations.  Distances  were  recorded 
between  the  center  of  a  tumor  in  these  images  and  in  the  phantom. 

D.  Localization 

Localization  was  assessed  in  a  forced  choice  task,  where  each  3.0-cm-diameter 
spherical  search  volume  Q  contained  one  and  only  one  tumor.  In  the  phantom,  the  tumor 
185  was  centered  in  Q  at  coordinates  (x,,  yt,  z,).  A  numerical  observer  sampled  Q  every  0.3 
cm  in  Cartesian  space.  The  voxel  with  the  greatest  response  and  its  26  neighbors  were 
then  sampled  more  finely  at  every  0. 1  cm.  The  observer  selected  subvoxel  location 
i(xm  ,ym  ,zm)  within  L>  of  noisy  image  g  that  produced  the  greatest  scalar  response  from 


un-normalized  cross  correlation  (XC)  with  an  observer  template  w: 


(2) 


i(xm,ym,zm)  =  arg  max  ( XCf ) , 

ieQ 

where 

xq  =  wgi.  (3) 

190  Thus,  location  (xm,  ym,  zm)  was  recorded  as  the  measured  tumor  position.  The  3D 

template  w  was  a  sphere  -  the  background-subtracted  tumor  signal  as  defined  in  the 
XCAT  activity  phantom.  Template  dimensions  were  1.5x1. 5x1. 5  cm3.  The  same 
template  w  was  used  for  all  12  tumor  sites  and  for  both  reconstruction  methods  and 
variations  of  cross  correlation.  Note  that  the  spherical  template  was  not  matched  exactly 
195  to  a  tumor  because  tumor  appearance  in  the  SPECT  images  was  distorted  due  to  the 

effects  of  distance-dependent  spatial  resolution,  iteration  number,  and  post-reconstruction 
smoothing. 

Normalized  cross  correlation  (NXC)  was  also  investigated  for  estimating  tumor  positions: 


NXC;-  = 


W’gj 


(4) 


Localization  error  dr  was  calculated  as  the  Euclidean  distance  between  a  true 
200  tumor  position  and  a  measured  tumor  position: 

dr  -  ■  ( 5 ) 

Localization  error  dr  is  always  positive  and  also  discrete  due  to  the  0.1 -cm  sampling  on  a 
Cartesian  grid.  Localization  errors  were  recorded  across  the  noisy  image  ensembles  for 
each  tumor  site.  Localization  errors  were  compared  for  DRC-XC,  DRC-NXC,  nDRC- 
XC,  and  nDRC-NXC  using  site-specific  optimal  smoothing  and  iteration  number 
205  parameters.  Differences  in  ensemble  localization  errors  were  assessed  for  statistical 
significance  using  the  Wilcoxon  signed-rank  test. 


For  each  tumor  site,  direction-dependent  localization  errors  were  computed  every 
3°  over  180°.  These  errors  were  indexed  by  an  angle  0.  Angle  #was  defined  in  a 
transaxial  plane  as  follows:  <9  was  equal  to  0  for  the  view  from  left  lateral  to  center  of 
210  rotation,  n/2  for  the  anterior  view,  and  7t  for  the  right  lateral  view.  Direction-dependent 
localization  bias  ( b ■*■&)  and  precision  indexed  to  a  specific  value  of  6,  were 

calculated  along  the  dimension  perpendicular  to  0.  As  such,  localization  precision  at  the 
anterior  view  described  the  lateral  spread  of  localization  errors,  and  localization  bias  -  the 
lateral  bias.  Figure  4  illustrates  this  convention  using  data  from  tumor  G.  The  direction- 
215  dependent  localization  errors  were  defined  mathematically  as  follows:  Measured  tumor 

coordinates  [xm  k,  ymk )  from  noisy  image  A:  and  true  tumor  coordinates  [xt,  yt )  were 
rotated  by  angle  fusing  matrix  T: 

cos 6  -sin# 

T  = 

sin  6  cos  6 


to  yield 


and 


(7) 


(8) 


220  For  the  direction  perpendicular  to  6,  direction-dependent  localization  bias  was 


calculated  as 


N=SO 
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and  direction-dependent  localization  precision  as 
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(ii) 


For  each  angle  6 direction-dependent  localization  precision  and  bias  were  analyzed  with 
225  respect  to  the  tumor-collimator  distance  and  also  with  reference  to  -  the  fraction 
of  photons  from  a  tumor  emitted  towards  the  detector  not  attenuated  by  the  phantom. 
These  metrics  are  related  to  distance-dependent  spatial  resolution  and  attenuation  of 
tumor  signal. 

III.  RESULTS 

230  A.  Reconstructed  Images 

Noisy  OSEM  images  were  smoothed  with  3D  Gaussian  kernels,  ranging  from  0  to 
2.5  cm  in  FWHM,  at  iterations  1-25.  Representative  nDRC  and  DRC  images, 
reconstructed  from  the  same  projection  set,  are  displayed  on  the  same  gray  scale  in 
Figure  5.  Anterior,  superficial  tumors  appear  hotter  in  DRC  images  than  in  nDRC 
235  images.  With  greater  smoothing,  the  heart  begins  to  appear  near  tumor  site  I.  Noise 
texture  differed  between  nDRC  and  DRC  images.  Figure  6  shows  ensemble-average 


images. 


Distances  were  recorded  between  tumor  positions  in  the  XCAT  phantom  and  in 
noise-free  background-subtracted  reconstructed  images.  For  DRC  images,  the  mean 
240  distance  and  standard  deviation  were  0.106  +/-  0.068  mm,  as  calculated  across  the  12 
tumor  sites  and  25  iterations.  The  mean  distance  was  0.128  +/-  0.060  mm  in  nDRC 
images.  These  shifts  in  tumor  positions  are  relatively  small  compared  with  the 
localization  errors  reported  in  this  study  below. 

B.  Localization 

245  The  mean  localization  error  was  calculated  from  the  80-image  ensemble  for  each 

tumor  site,  reconstruction  method,  variation  of  cross  correlation,  iteration  number,  and 
degree  of  smoothing.  Mean  localization  errors  are  mapped  out  as  a  function  of  these 
parameters  in  Figure  7  where  darker  color  signifies  better  localization.  A  white  crossed 
box  marks  the  combination  of  iteration  number  and  smoothing  that  resulted  in  the  lowest 
250  mean  localization.  Subsequent  comparisons  were  performed  using  these  optimal 

parameter  values.  In  Table  I,  optimal  combinations  of  smoothing  and  iteration  number 
are  listed  with  corresponding  mean  localization  errors. 

Localization  errors  -  calculated  across  the  80-image  ensemble  using  Eq.  (5)  -  are 
summarized  in  Figure  8.  The  four  combinations  of  reconstruction  method  and  cross 
255  correlation  -  nDRC-XC,  nDRC-NXC,  DRC-XC,  DRC-NXC  -  are  represented  by  four 
boxplots  at  each  tumor  site.  Each  box  plot  includes  80  localization  estimates.  Also 
shown  is  a  box  plot  representing  random  chance  -  the  uniform  probability  of  selecting 
any  discrete  0.1-cm-spaced  location  within  a  3-cm-diameter  search  volume  centered  on 
tumor.  The  random  chance  boxplot  shows  that  more  than  75%  of  potential  tumor 
260  locations  are  9  mm  or  farther  from  a  true  tumor  position  and  thereby  serves  as  a  useful 


benchmark  for  gauging  the  amount  of  information  provided  by  SPECT  for  target 
localization.  Note  that  localization  appears  worse  than  expected  due  to  random  chance  at 
tumor  site  I  using  un-normalized  cross  correlation  because  the  numerical  observer  selects 
the  much  hotter  heart. 

265  Ensemble  localization  errors  are  described  more  completely  by  plotting  the 

fraction  of  correct  localizations  within  a  radial  distance  from  a  true  tumor  position,  as 
illustrated  in  Figure  9.  Two  tumor  sites  -  one  each  of  relatively  good  and  relatively  poor 
localization  -  are  plotted  together.  For  reference,  and  as  described  above,  the  random 
chance  curve  is  plotted  in  black.  In  these  plots,  localization  is  better  for  a  curve  that  is 
270  above  and  to  the  left  of  another. 

The  Wilcoxon  signed-rank  test  was  used  to  assess  for  significant  differences  in 
ensemble  localization  errors  at  each  tumor  site  between  nDRC-XC,  nDRC-NXC,  DRC- 
XC,  and  DRC-NXC.  P-values  are  reported  in  Table  II,  with  values  less  than  0.05  in  bold. 
Localization  was  the  worst,  and  significantly  so,  using  nDRC-XC  at  tumor  sites  A,  B,  C, 
275  D,  and  F.  At  these  sites,  mean  localization  was  best  using  DRC-NXC  (Table  I);  however, 
localization  differences  were  not  statistically  significantly  different  than  DRC-XC  and 
nDRC-NXC  at  each  site.  Near  the  heart  at  site  I,  NXC  was  significantly  better  than  XC. 
Localization  performance  was  best  using  nDRC  for  sites  H,  J,  and  L. 

In  Figure  10,  tumor  position  estimates  are  overlaid  on  the  attenuation  map.  There 
280  are  80  markers  at  each  tumor  site  encompassing  estimates  from  the  image  ensemble. 

Localization  error  patterns  depend  largely  on  the  tumor  site.  For  instance  at  site  B  in  the 
right  breast,  position  estimates  are  tightly  clustered  near  the  true  tumor  position.  In 
contrast,  position  estimates  are  distributed  more  diffusely  at  site  G  in  the  right  posterior 


rib  where  there  was  a  noticeable  anisotropy.  Note  the  similarity  between  localization 
285  patterns  in  Figure  10  and  the  appearance  of  reconstructed  tumors  in  Figure  6. 

Localization  was  performed  in  3D  images  reconstructed  from  all  detector  views. 
The  quality  of  information  from  specific  detector  views  was  predictive  of  direction- 
dependent  localization  error.  This  effect  is  illustrated  in  Figure  1 1  where  position 
estimates  at  site  G  are  plotted  with  reference  to  the  detector  trajectory.  Target 
290  localization  was  more  precise  when  the  detector  was  in  close  proximity  to  tumor  as 
compared  with  a  more  distal,  orthogonal  view. 

Anisotropies  in  localization  patterns  were  quantitatively  analyzed  by  calculating 
direction-dependent  localization  bias  b±d  and  precision  as  a  function  of  angle 
Q,  using  Eqs.  (9,10).  Directional  dependencies  of  b±d  and  &-&  are  summarized  for  each 
295  tumor  site  in  Figure  12.  Beneath  localization  bias  and  precision  plots  are  two  additional 
curves:  One  curve  marks  the  distance  between  tumor  position  and  patient-side 
surface  of  the  collimator  at  6.  The  other,  filled  curve  shows  attenuation  survival 
probability  y1#  along  6.  Attenuation  survival  probability  is  the  fraction  of  photons 
emitted  from  the  tumor  towards  the  detector  not  attenuated  by  the  phantom.  Direction- 
300  dependent  localization  bias  and  precision  tend  to  be  best  along  directions  0  with  low 
values  of  d^®  and  high  values  of  particularly  at  sites  A-F.  Values  of  d^®  and 
y^were  averaged  over  the  detector  trajectory  and  are  reported  in  Table  III. 

IV.  DISCUSSION 

In  a  previous  computer-simulation  study,  we  investigated  localization 
305  performance  in  SPECT  images  as  a  function  of  scan  time,  tumor  size,  tumor-to- 
background  uptake  ratio,  and  tumor  position.2  In  the  present  study,  localization 


performance  was  compared  for  two  reconstruction  methods  -  nDRC  and  DRC  -  and  two 
variations  of  cross  correlation  -  XC  and  NXC  -  between  the  observer  template  and 
SPECT  images.  In  addition,  localization  error  patterns  were  analyzed  for  directional 
310  dependencies.  Metrics  were  developed  to  describe  direction-dependent  localization  bias 
and  precision.  These  directional  errors  were  related  to  the  detector  trajectory. 

In  this  study,  localization  performance  was  assessed  using  numerical  observers 
that  had  no  knowledge  of  noise  correlations.  As  such,  the  observers  were  suboptimal. 
Prewhitening  of  noise  correlations  may  improve  localization  performance.  Though 
315  prewhitening  presents  practical  challenges,  methods  for  rapidly  computing  the  noise 

covariance  matrix  and  for  using  channels  to  reduce  covariance  dimensions  are  areas  of 
active  research.17'19  That  noted,  for  the  detection  of  hot  spheres  in  unknown  locations, 
non-channelized,  non-prewhitening  matched  filters  have  been  shown  to  outperform 
human,  channelized  non-prewhitening,  and  channelized  Hotelling  observers.  Though 
320  channelized  observers  can  better  match  human  performance  by  using  only  information  in 
certain  frequency  channels  20 ,  it  was  not  our  objective  here  to  match  human  performance. 

The  non-prewhitening  matched  filter  can  be  defined  as  wnpw  =  gi  -  go,  where  g, 
and  g0  are  expected  values  for  tumor-present  and  tumor-absent  images  and  wNPW  -  the 
expected  difference  -  is  an  image  of  the  tumor.21  Thus,  the  non-prewhitening  matched 
325  filter  accounts  for  location-dependent  distortions  of  the  tumor  appearance  in  a  SPECT 

image.  In  practice  clinically,  the  matched  filter  may  not  be  obtainable.  In  recognition  of 
this  limitation,  the  template  used  in  this  study  approximated  the  matched  filter  with  a 
sphere.  Localization  performance  might  be  different  using  the  non-prewhitening 
matched  filter  or  another  observer  template  instead  of  a  spherical  template.  In  some 


330  clinical  scenarios,  it  may  be  that  a  tumor  template  which  better  approximates  the  matched 
filter  could  be  estimated  from  a  treatment  planning  SPECT  image. 

Comparisons  in  localization  performance  were  made  between  nDRC-XC,  nDRC- 
NXC,  DRC-XC,  and  DRC-NXC  for  the  iteration  number  and  degree  of  smoothing  that 
yielded  the  lowest  mean  localization  error.  Post-reconstruction  smoothing  of  SPECT 
335  images  with  Gaussian  kernels  has  been  shown  to  improve  LROC  performance.  In  this 
study,  smoothing  improved  localization  in  38  of  the  48  cases,  as  evaluated  across  the  12 
tumor  sites  and  the  4  combinations  of  reconstruction  method  and  variation  of  cross 
correlation.  At  each  tumor  site,  smoothing  was  always  helpful  for  NXC.  The  degree  of 
smoothing  optimal  for  localization  was  largely  determined  by  the  type  of  cross 
340  correlation  -  NXC  versus  XC.  NXC  localization  was  optimized  by  Gaussian  kernels  of 
approximately  three  times  greater  width  than  those  optimal  for  XC.  This  can  be 
understood  as  follows:  Smoothing  causes  activity  outside  of  the  search  volume,  or  lack 
thereof  for  tumors  near  body  contours,  to  be  blurred  into  the  search  volume  thus  creating 
activity  gradients.  XC  -  the  (un-normalized)  inner  product  -  is  biased  by  such  gradients, 
345  much  more  so  than  NXC.  For  example,  NXC  localization  is  much  better  than  XC 

localization  at  tumor  site  I,  which  is  near  the  heart.  By  allowing  for  greater  smoothing 
than  XC,  NXC  may  also  minimize  the  effects  of  high  frequency  noise.  Iteration  number 
was  optimized  in  conjunction  with  smoothing.  A  greater  number  of  image  updates 
(product  of  the  number  iterations  and  number  of  subsets)  were  needed  for  optimal 
350  localization  in  DRC  than  in  nDRC  images.  These  results  are  consistent  with  the  findings 
of  other  SPECT  studies.3, 23  In  two  cases,  localization  performance  was  best  at  iteration 
25.  In  the  absence  of  normalization,  the  XC  observer  selected  the  much  hotter  heart 


instead  of  tumor  I.  It  is  possible  that  at  higher  iteration  numbers,  cardiac  activity  would 
be  pulled  out  of  the  search  volume,  thereby  resulting  in  better  localization  performance. 
355  At  superficial  tumor  sites  A-F,  which  were  proximal  to  the  detector  trajectory, 

localization  performance  was  either  comparable  or  best  using  DRC-NXC.  However, 
differences  in  localization  results  were  not  statistically  significantly  different  from  DRC- 
XC  and  nDRC-NXC  at  each  tumor  site.  At  all  of  these  6  sites  except  for  E,  localization 
was  worst,  and  significantly  so,  using  nDRC-XC.  At  the  other  6  sites  where  attenuation 
360  was  2.9  times  more  severe  and  the  collimator  1.4  times  farther  from  tumor  on  average, 
DRC  typically  did  not  result  in  better  localization  than  nDRC.  Though  DRC  improved 
localization  at  G,  localization  was  worse  using  DRC  for  lung  tumors  H  and  J,  and  for  L 
near  the  spine.  These  differences  were  statistically  significant.  Study  results  support 
using  DRC  and/or  NXC  for  localization  of  proximal,  superficial  tumors.  This  study  did 
365  not  compare  DRC  and  nDRC  images  for  the  task  of  determining  the  size  or  extent  of 

tumor,  though  such  a  comparison  is  relevant  since  certain  biological  targets,  e.g.,  hypoxia, 
can  vary  with  time.24 

Tumor  proximity  to  the  detector  trajectory  had  a  major  effect  on  localization 
errors,  generally  more  so  than  the  reconstruction  method  or  type  of  cross  correlation. 

370  Close  detector  proximity  improves  spatial  resolution,  and  though  attenuation  only 

depends  on  the  material  between  the  detector  and  tumor,  it  too  is  typically  best  at  views 
of  close  detector  approach.  At  site  K,  mean  localization  error  was  very  poor;  however,  K 
was  at  a  comparable  depth  in  the  phantom  as  other  tumors,  e.g.,  E  and  G,  that  had  vastly 
better  localization  because  of  closer  proximity  to  the  detector  trajectory.  Although  the 


375  couch  will  limit  detector  approach  and  attenuate  photons,  localization  of  K  may  be 
improved  by  using  a  posterior  detector  trajectory. 

Localization  error  patterns  were  anisotropic.  Close  detector  proximity  to  a  tumor 
and  minimal  attenuation  were  good  predictors  of  a  direction  with  the  best  localization 
bias  and  precision.  The  direction-dependent  localization  pattern  at  site  G  is  interesting  in 
380  that  localization  is  relatively  good  in  one  dimension  but  very  poor  orthogonally.  This 
information  is  not  conveyed  when  only  using  distances  to  describe  localization  errors. 
Strong  directional  dependencies  have  important  implications  for  external  beam  radiation 
therapy  since  radiation  is  often  delivered  from  multiple  gantry  angles. 

The  optimal  iteration  numbers  and  smoothing  parameters  used  in  the  above 
385  comparisons  may  not  be  known  in  clinical  applications.  Figure  7  shows  that  there  are 
often  broad  regions  across  the  2D  iteration-number  /  smoothing-parameter  space  over 
which  localization  accuracy  is  near  to  its  optimal  value.  This  may  facilitate  successful 
application  of  practical  clinical  techniques  such  as  identifying  an  iteration  number  and 
degree  of  smoothing  that  work  well  across  a  broad  range  of  patients  and  anatomical 
390  locations,  perhaps  with  some  parameterization  based  on  factors  such  as  patient  diameter, 
tumor  depth,  and  tumor-to-background  uptake  ratio.  That  noted,  clinical  localization  may 
not  generally  be  able  to  achieve  the  optimal  performance  used  in  the  above  comparison 
studies.  We  note  that  other  applications  of  SPECT,  such  as  quantification  and  detection, 
are  also  often  performed  clinically  using  suboptimal  smoothing  and  numbers  of  iterations. 
395  This  study  used  an  ideal  attenuation  map.  If  implemented  clinically,  attenuation 

maps  could  be  derived  from  an  offline  CT,  registered  with  SPECT  or  CBCT,  or  from  on¬ 
board  MV  or  kV  CBCT.  Because  these  x-ray  beams  are  polychromatic  and  "mTc  is 


monochromatic,  the  SPECT  attenuation  map  cannot  be  measured  directly,  but  it  can  be 
calculated.  CBCT  images  are,  however,  subject  to  artifacts  caused  by  insufficient 
400  sampling,  beam  hardening,  and  scatter.  Inaccuracies  in  attenuation  maps  can  cause  errors 
in  SPECT  images.26  Further  investigation  is  needed  to  determine  the  effect  of  an 
imperfect  attenuation  map  on  localization  performance.  The  simulated  detector  had  ideal 
scatter  rejection.  Scatter  has  been  shown  to  degrade  LROC  performance,  though  this 

97 

effect  can  be  lessened  by  compensating  for  scatter  during  image  reconstruction. 

405  V.  CONCLUSIONS 

This  computer-simulation  study  compared  localization  performance  for 
normalized  and  un-normalized  numerical  observers,  which  were  used  to  estimate  tumor 
positions  in  SPECT  images,  reconstructed  with  and  without  DRC.  For  tumors  localized 
to  <  2  mm  on  average,  which  are  good  candidates  for  SPECT-guided  radiation  therapy, 
410  localization  performance  typically  improved  by  compensating  for  the  detector  response 
function  and  by  using  a  normalized  observer.  Localization  errors  were  analyzed  across 
the  image  ensemble  and  were  found  to  be  anisotropic  and  dependent  on  tumor  position 
relative  to  the  detector  trajectory.  Directional  dependencies  are  attributable  to  the  effects 
of  attenuation  and  distance-dependent  spatial  resolution.  Direction-dependent 
415  localization  error  patterns  have  important  implications  for  radiation  therapy  and  are 
relevant  to  SPECT  imaging  in  general. 
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FIGURE  CAPTIONS 


425  Fig.  1.  (a)  Simulated  activity  ratios  were  20:6:1  in  heart,  tumor,  and  other  tissues,  (b) 
Twelve  tumor  sites,  labeled  A-L,  were  in  axial  slice  superior  to  heart. 

Fig.  2.  (a)  LEHR  collimator  hole  geometry  -  length  l  of  2.635  cm  and  diameter  d  of 
0. 140  cm  -  determines  angles  <p  through  which  rays  may  pass  directly  to  the  detector,  (b) 
430  Cones  of  49  rays  were  used  to  simulate  distance-dependent  collimator  spatial  resolution, 
(c)  Shows  agreement  between  simulated  and  measured  line-source  profiles  at  distances  of 
15  to  30  cm.  Measurements  were  made  with  a  Trionix  Triad  scanner  using  LEHR 
collimation. 

435  Fig.  3.  The  conformal  180°  detector  trajectory  is  overlaid  on  the  attenuation  map.  At  the 
60  detector  views,  a  circle  marks  the  collimator  position.  The  detector  width  is  indicated 
by  gray  bars  at  the  right  lateral  and  anterior  views  where  distances  to  the  axis  of  rotation 
are  25.7  cm  and  23.9  cm  respectively. 

440  Fig.  4.  (a)  An  anisotropic  localization  error  pattern  is  marked  by  filled  diamonds,  (b) 
Direction-dependent  localization  precision  &-&  is  indexed  by  6.  Note  that  6  is  also 
shown  above. 

Fig.  5.  Representative  images  (nDRC,  DRC)  reconstructed  using  the  same  projection  data 
445  are  displayed  as  function  of  iteration  number  and  degree  of  post-reconstruction 
smoothing. 

Fig.  6.  (a)  nDRC  and  (b)  DRC  images  were  averaged  over  the  80-image  ensemble  at 
iteration  number  10.  Spherical  tumors  appear  elliptical  and  vary  in  orientation  and 
450  degree  of  eccentricity.  Distortions  are  related  to  tumor  position  relative  to  the  detector 
trajectory. 

Fig.  7.  Mean  localization  error  dr  for  tumor  sites  A-L  is  plotted  as  a  function  of 
smoothing,  iteration  number,  reconstruction  method,  and  cross  correlation.  Within  each 
455  broad  column  and  row,  mean  localization  error  is  displayed  as  a  function  of  iteration 

number  (1-25)  on  the  horizontal  axis  and  smoothing  (0  to  2.5-cm-FWHM)  on  the  vertical 
axis.  Labeled  axes  are  shown  above  the  phantom,  using  tumor  F,  DRC-XC  as  an 
example.  The  combination  of  iteration  number  and  smoothing  with  the  best  mean 
localization  error  is  marked  by  a  white  crossed  box. 

460 

Fig.  8.  Box  plots  summarize  the  ensemble  localization  errors  of  tumors  A-L.  Black 
horizontal  bars  denote  medians,  filled  rectangles  the  interquartile  range  (IQR),  whiskers 
the  last  data  point  within  1 .5  x  IQR  from  the  25th  and  75th  percentiles,  and  filled  markers 
the  outliers.  Marked  by  *  is  a  boxplot  representing  random  chance. 

465 

Fig.  9.  The  fraction  of  correct  localizations  is  plotted  as  a  function  of  radius  from  the  true 
tumor  location.  Dashed  lines  represent  normalized  cross  correlation,  solid  lines  un- 


normalized  cross  correlation.  DRC  lines  are  lighter  than  nDRC  counterparts. 

Localization  by  random  chance  is  shown  by  a  black  curve. 

470 

Fig.  10.  Estimates  of  tumor  positions  from  the  image  ensembles  are  plotted  on  the 
attenuation  map.  For  reference,  tumor  diameters  are  marked  by  solid  white  lines  and 
search  region  boundaries  by  dotted  white  lines.  Correct  localizations  are  centered  exactly 
in  these  circles. 

475 

Fig.  1 1 .  (a)  The  detector  is  highlighted  at  two  views  -  one  proximal  to  tumor  and  the 
other  more  distal  and  orthogonal,  (b)  Position  estimates  at  site  G  show  that  more  precise 
localization  is  associated  with  a  more  proximal  detector  view  where  spatial  resolution 
was  relatively  good  and  attenuation  less  severe. 

480 

Fig.  12.  (a)  Direction-dependent  localization  bias  and  precision  are  plotted  as  a  function 
of  <9  for  tumors  A-F.  Horizontal  dashed  lines  mark  ideal  bias  and  precision.  Also  plotted 
are  curves  showing  tumor-collimator  proximity  and  the  survival  probability  for  photons 
from  a  tumor  to  escape  the  phantom  without  being  attenuated,  (b)  Shows  results  for 
485  tumors  G-L.  Direction-dependent  localization  bias  curve  is  truncated  at  site  I. 

TABLE  CAPTIONS 

Table  I.  Mean  localization  errors  are  reported  with  optimal  combinations  of  smoothing 
490  and  iteration  number. 

Table  II.  Differences  in  ensemble  localization  errors  were  assessed  using  Wilcoxon 
signed-rank  test.  Bold  indicates  p  <  0.05. 

495  Table  III.  Average  tumor-collimator  distance  d±d  and  attenuation  survival  probability  f~0 
are  reported  for  tumors  A-L. 
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Quantitative  Analysis  of  a  Compact  On-board  SPECT  Detector  for  ROI  Reconstruction 

JR  Roper,  JE  Bowsher,  FF  Yin 

Purpose:  Single  photon  emission  computed  tomography  (SPECT)  imaging  on-board  radiation  therapy  machines  may 
enhance  radiation  therapy.  Use  of  a  light-weight,  compact  detector  would  be  desirable  because  of  maneuverability,  yet 
is  challenged  by  truncation  and  consequently  reconstruction  artifacts.  We  hypothesize  that,  for  a  sizeable  volume 
surrounding  the  tumor  target,  and  with  properly  chosen  detector  trajectories,  compact  detectors  can  provide  image 
quality  that  is  comparable  to  that  of  full-size  SPECT  detectors. 

Methods  and  Materials:  On-board  SPECT  imaging  was  computer  simulated  for  detectors  with  active  surface  widths 
of  21 .2  and  40.0  cm.  Detector  trajectories  were  selected  such  that  the  common  volume  encompassed  an  8-cm- 
diameter  ROI  surrounding  each  tumor  in  the  torso  of  a  NCAT  phantom.  For  deep  tumors,  two  common  volumes  - 
interior  and  partially  exterior  -  were  investigated.  Radiotracer  distribution  was  modeled  for  99mTc-Sestamibi.  Noise- 
free  and  noisy  projection  images  were  generated  with  an  analytical  simulator  that  models  non-uniform  attenuation, 
collimator  &  detector  efficiency,  and  spatially-varying  spatial  resolution.  Images  were  reconstructed  by  OSEM.  In 
reconstructed  images,  root  mean  square  error  and  recovered  tumor  activity  were  analyzed  as  a  function  of  detector 
width  and  detector  trajectory. 

Results:  Image  truncation  reduced  overall  image  quality.  For  example,  RMS  error  over  an  entire  noise-free  image 
slice  was  18  -  35%  worse  with  the  21-cm-wide  detector  versus  the  larger  detector.  However,  for  regions  within  the 
common  volume  for  both  detector  widths,  RMS  errors  differed  by  less  that  2%.  Similar  results  were  observed  for  noisy 
images.  Moreover,  the  recovered  fraction  of  tumor  activity  was  comparable,  except  for  when  the  smaller-detector 
common  volume  was  interior  -  the  recovered  activity  was  reduced  by  1 1  %. 

Conclusions:  Preliminary  results  show  that  an  8-cm-diameter  ROI  can  be  reconstructed  using  21  or  40-cm-wide 
detectors  with  comparable  RMS  errors  and  recovered  tumor  activities,  supporting  the  proposal  that  on-board  target 
localization  could  be  accomplished  using  compact  SPECT  detectors. 
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Are  post-RT  cardiac  perfusion  defects  due  to  cardiac  toxicity, 
or  are  they  artifacts  from  attenuation  changes  in  surrounding  soft  tissues? 

JR  Roper,  AA  Manzoor,  JB  Bowsher,  FF  Yin,  S-M  Zhou,  TZ  Wong,  S  Borges-Neto,  JL  Hubbs,  S  Demirci,  LB  Marks. 


Purpose:  Following  RT  for  left-sided  breast  cancer,  perfusion  defects  (i.e.  cold  regions)  are  often  noted  on  SPECT  images. 
Almost  all  new  defects  are  in  the  cardiac  apex.  Since  SPECT  is  not  absolutely  quantitative,  defects  are  detected  by 
comparing  regional  intensities  on  pre-RT  and  post-RT  images.  Further,  as  SPECT  images  are  not  corrected  for  attenuation, 
and  since  RT  may  increase  soft  tissue  density,  apparent  defects  may  be,  at  least  partially,  attributable  to  changes  in 
attenuation  of  surrounding  soft  tissue.  We  herein  perform  a  series  of  computer-based  simulations  to  assess  this  possibility. 

Methods:  The  perfusion  tracer  "mTc-Sestarnibi  was  simulated  in  a  female  NCAT  phantom.  Two  tangential  breast  fields 
were  considered:  the  deep  border  of  one  excluded  the  heart  while  a  second  included  2.5  cm  of  anterior  heart.  Within 
RT  fields,  tissue  densities  were  increased  in  increments  of  10%  up  to  50%.  Noise-free  images  were  simulated  for  a 
cardiac  SPECT  protocol  -  180°  circular  orbit,  60  views,  parallel  hole  collimation  -  by  modeling  photon  detection 
efficiency,  distance-dependent  spatial  resolution,  and  non-uniform  attenuation.  Images  were  reconstructed  without 
attenuation  correction  by  OSEM  (10  subsets,  20  iterations)  onto  a  grid  with  0.36  cm  wide  voxels.  For  varying  degrees 
of  change  in  soft  tissue  density,  apparent  reductions  within  the  apex  were  noted,  and  compared  with  effects 
elsewhere.  Simulations  were  repeated  with  the  left  (222  ml)  breast  removed  to  model  mastectomy. 


Results:  As  tissue  density  is  increased  by  10,  20,  30,  40,  and  50%,  corresponding  reductions  in  apparent  apical  perfusion, 
for  the  different  scenarios,  are,  respectively: 

13,  24,34,  42,  50%; 

8,  15,  22,  28,  34%; 

10,  19,  26,  34,40%; 

5,  9,14,18,21%. 


intact  breast,  2.5  cm  heart  in  field 
intact  breast,  heart  sparing 
mastectomy,  2.5  cm  heart  in  field 
mastectomy,  heart  sparing 


Changes  in  cardiac  perfusion  were  not  uniform;  the  apex  was  affected  more  than  other  areas. 


Conclusion:  RT-induced  increases  in  soft-tissue  attenuation  can  cause  apparent  changes  in  apical  perfusion  that  might 
account  for  some  of  the  defects  seen  clinically.  However,  the  simulations  also  suggest  that  apparent  defects  should  be  seen 
even  if  the  heart  is  excluded  from  the  RT  beam.  Since  perfusion  defects  are  rarely  seen  in  patients  with  no  heart  in  the  field, 
the  possible  changes  in  soft  tissue  attenuation/density  are  likely  not  the  primary  cause  of  perfusion  defects  seen  clinically. 
Additional  work  is  needed  to  understand  the  degree  and  impact  of  RT-induced  soft  tissue  changes. 


Appendix  F 


On-board  SPECT  Imaging:  Improved  Target  Localization  via  Spatial  Resolution  Compensation 

JR  Roper,  JE  Bowsher,  FF  Yin 

Purpose:  Single  photon  emission  computed  tomography  (SPECT)  is  being  investigated  for  imaging  on-board  radiation 
therapy  machines  to  localize  biological  targets.  SPECT  spatial  resolution  degrades  as  distance  increases  between  the 
detector  and  radiotracer,  an  effect  of  collimation  which  results  in  spatially-varying,  asymmetric  blurring.  Here  we  assess 
localization  accuracy  in  noisy  images  that  were  reconstructed  with  and  without  spatial  resolution  compensation. 

Methods:  Uptake  of  "mTc-Sestamibi  was  simulated  in  a  female  NCAT  phantom.  Twenty,  1 .4-cm-diameter  biological 
targets  were  simulated  in  the  torso.  Target  activity  ratios  were  3:1  and  6:1  relative  to  a  background  concentration  of 
0.25  pCi/ml.  Spatially-varying  spatial  resolution  was  modeled  in  projection  images  for  a  half-circular  anterior  orbit.  An 
ensemble  of  25  noisy  projection  images  was  degraded  with  Poisson  noise  to  simulate  8  minute  scans.  Images  were 
reconstructed  using  10-subset,  5-iteration  OSEM  with  and  without  spatial  resolution  compensation.  Compensation 
involved  tracing  cones  of  49  rays  through  the  collimator  acceptance  angle,  as  opposed  to  tracing  1  central  ray. 
Nonprewhitening  matched  filters  were  used  to  measure  target  centroids  within  3D,  2.9-cm-diameter  search  volumes. 
Localization  accuracy  was  calculated  as  the  mean  distance  between  measured  and  true  target  centroids  over  the 
ensemble. 

Results:  Localization  accuracy  was  less  than  4  mm  for  6  targets  using  either  reconstruction  method.  For  these  targets, 
spatial  resolution  compensation  always  improved  localization  accuracy.  On  average,  localization  accuracy  improved  from  2.5 
mm  to  1 .9  mm.  The  maximum  improvement  was  1  mm.  With  spatial  resolution  compensation,  two  targets  were  localized  with 
1  mm  accuracy. 

Conclusion:  Preliminary  results  show  improved  localization  in  on-board  SPECT  images  with  spatial  resolution  compensation 
during  reconstruction.  These  findings  are  important  for  assessing  the  feasibility  of  on-board  SPECT  imaging  and  warrant 
further  investigation. 


Appendix  G 


On-board  SPECT:  Localization  Study  using  Scanner-acquired  Data 

JR  Roper,  JM  Wilson,  JE  Bowsher,  TG  Turkington,  FF  Yin 

Purpose:  In  previous  computer-simulation  studies  we  have  shown  that  on-board  SPECT  imaging  is  a  promising 
approach  to  localizing  biological  structure.  The  purpose  of  this  study  is  to  investigate  SPECT  localization  using 
scanner-acquired  data. 

Method  and  Materials:  A  torso  phantom  was  filled  with  0.75pCi/ml  of  "mTc.  Five  targets  (A-E)  of  6:1  activity  ratio  - 
diameters  of  22  to  34mm  -  were  positioned  in  phantom.  SPECT  data  were  acquired  with  one  Trionix  Triad  detector. 
The  product  of  activity  and  1.5-minute  count  time  models  a  4.5-minute  clinical  scan.  Ensembles  of  10  images  were 
reconstructed  using  MLEM  for  3  detector  trajectories:  360°,  180°  left  lateral,  and  180°  right  lateral.  Localization  was 
assessed  across  the  noisy  ensembles  using  non-prewhitening  observers  that  were  forced  to  select  the  most 
suspicious  target  location  from  a  3.6-cm  diameter  search  region  centered  on  target.  True  centroids  were  derived  from 
a  12-hour  SPECT  scan  that  registered  with  CT  to  within  1/3mm.  Localization  error  was  calculated  as  the  distance 
between  true  and  measured  centroids.  Significance  was  assessed  with  the  Wilcoxon  rank  sum  test. 

Results:  Mean  localization  errors  for  targets  (A-E)  are  A)  0.80mm,  B)  1 ,2mm,  C)  1 ,9mm,  D)  1 ,2mm,  and  E)  0.90mm 
using  the  180°  trajectory  more  proximal  to  each  target.  With  the  360°  trajectory,  mean  errors  are  A)  1 .3mm,  B)  1 .3mm, 
C)  2.1mm,  D)  2.2mm,  and  E)  1.5mm.  Differences  between  the  above  180°  and  360°  trajectories  are  significant  (p  < 
0.05)  at  A,  D,  and  E.  Mean  error  is  typically  much  worse  using  the  more  distal  180°  trajectory:  A)  4.5mm,  B)  3.1mm,  C) 
2.8mm,  D)  4.4mm,  and  E)  6.9mm. 

Conclusion:  This  study,  using  scanner-acquired  SPECT  data,  demonstrates  mean  localizations  to  within  1  or  2mm  with 
4.5-minute  scans.  The  study  also  shows  that  detector  trajectory  can  substantially  affect  target  localization  and  should 
be  optimized  for  specific  target  sites. 


